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INTRODUCTION 


The purpose of the computer program GRADE is to numerically solve the 
differential equations that describe heat pipes with graded-porosi ty 
fibrous wicks. Such wicks have an axial variation in porosity so that, 
at the maximum heat- transfer rate, the porosity is just low enough for 
the wick to support the local liquid-flow pressure drop plus any hydro- 
static head. Thus, the highest possible permeability is obtained along 
the wick. For comparison of graded- to uniform-porosity wicks, GRADE can 
also compute the performance of the latter. In fact, the user can have 
GRADE calculate the porosity for a given fiber diameter, or both the poros- 
ity and fiber diameter, that gives the highest heat-transfer capacity for 
a uni form-porosi ty wick. 

In the graded-porosi ty case, the user specifies an initial porosity* 
at the condenser end and either a final porosity or maximum liquid stress** 
at the evaporator end. GRADE then calculates the optimum porosity varia- 
tion along the wick and the maximum heat-transfer rate. A summary of 
GRADE'S capabilities follow: 

o Calculation of optimum porosity variation and corresponding 
maximum heat- transfer rate. 

o Calculation of maximum heat-transfer rate at other than the 
wick's design condition. 

o Both earth - gravity and zero-gravity computations. 

o Heat pipes having multiple sections each having different tilts. 

o Multiple heat input and output zones. 


*If the user wishes, he can have GRADE compute the initial porosity as the 
highest for which the wick's condenser end will self-prime in a gravita- 
tional field. 

**We define stress as the local vapor-liquid pressure difference. 
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0 Calculation of the amount of liquid contained in the wick. 

o Calculation of the performance of uniform-porosity wicks with 
a specified constant porosity. 

o Calculation of optimum uni form-porosity wicks. 

The mathematical model for flow through a wick is based on the TRW 
Report No. 21160-6005-RU-00 , "Capillary Flow through Partially Saturated 
Wicks", however that analysis has been corrected and modified in light of 
recent theoretical and experimental research. A brief discussion of the 
mathematical basis for GRADE follows. 

CAPILLARY PRESSURE 

The model of capillary pressure generated by a fibrous wick includes 
hysteresis, that is, the capillary pressure generated for filling being 
significantly lower than that for emptying. The expression for the capil- 
lary pressure P generated by a wick of porosity <p and of diameter 6 
is given by 

P = 3.247RH LbiL _£_ . 0) 

<P 0 

Here, o is the surface tension and H is an empirical constant equal to 
unity if the wick is filling and 1.955 if emptying. The empirical constant 
R is the ratio of the maximum stress of a full region to the larger stress 
necessary to half empty the region. Its value is approximately 0.6 . 

Eq. (1) is derived by first calculating the component of force perpen- 
dicular to a vapor-liquid interface exerted by surface tension acting on a 
single fiber that penetrates the interface at an arbitrary orientation. 

The force components of all fibers penetrating a given area are then summed 
with proper account taken of the statistical distribution of their orient- 
ations, and the sum is equated to the product of the capi 11 ary pressure 
and the liquid area. 

Eq. (1) can be solved for <j> , which results in 

3 . 247RHo/6 

* P + 3.247RHa/6 (2) 
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PERMEABILITY 


The permeability K is the proportionality factor in Darcy's law, 

, dP 

which relates the superficial mass flux M to the pressure gradient 
and the kinematic viscosity u : 


The expression we use for the permeability in terms of the porosity $ and 
the fiber diameter 5 is 


K = - d — 5L 


ii 

* |40-*)-(l-*) 2 -2 ln(W)-3 


-1 


ln{l-$) 


HM ) 2 

1 +( 1 -<J>) ^ 


(3) 


The empirical constant D is introduced to obtain close agreement between 
theory and experiment. 

This expression is derived by first calculating the drag on an arbi- 
trary orientated fiber segment in the presence of neighboring fibers. 

The drags on all such segments contained in a cylindrical region whose 
axis is aligned with the flow direction are summed with account taken of 
the statistical distribution of their orientations. The total drag is 
then equated with the difference in pressure forces at each end of the 
cylinder. 


OPTIMUM UNIFORM WICK 

With the use of Eq. (1) for the capillary pressure and Eq. (3) for 
the permeabi 1 i ty , we can optimize the capacity of a uniform-porosi ty wick 
with respect to its fiber diameter and porosity. The results for operation 
in a gravitational field are 

o The optimum porosity is 0.79 and the optimum fiber diameter 
is such that Eq. (1) gives twice the hydrostatic contribution 
to the stress at the evaporator end. 

o If the fiber diameter is fixed, the optimum porosity is given 
by Eq. (2) with P equal to twice the hydrostatic contribution 
to the stress at the evaporator end. 
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MATHEMATICAL PROBLEM STATEMENT 

In this formulation, we neglect vapor-flow pressure drop. Three 
first-order differential equations govern the stress P (at the top of 
the wick if in a gravi tational field), the mass flow rate m , and the 
cumulative mass M of liquid in the : wick: 



( 4 ) 

( 5 ) 


= pA 4>(x ) , W 

J / « 

where u,p,g,x,A,g^- 3 and Q(x) are, respectively, the 
kinematic viscosity, density, gravitational acceleration, latent heat, 
wick area, heat-pipe slope, and heat input per unit length. The perme- 
ability K is given by Eq. (3), and for the optimum graded-porosity wick, 
the porosity <|>(x) is set just low enough to generate the required capil- 
lary pressure. In particular, <j>(x) is taken as the lowest of: 


(1) The porosity required to sustain the stress P , which is 
calculated from Eq. (2) with the hysteresis constant H = 
1.955 . 

(2) The initial porosity , which, if not specified, is 
calculated from Eq. (2) with P = pgh w ( h w is the wick 
height) and H = 1 . 

(3) The local self-priming porosity, which is calculated from 
Eq. (2) with H = 1 and P = pg[h(x) + h w ] . 

In zero gravity, priming is not a factor, and the porosity <|>(x) is 
the lower of Items (1) and (2) above. In non-zero gravity, if the user does 
not require the wick to self-prime along its entire length at the operating 
elevation, then the self-priming porosity given by Item (3) is not 
considered. For this case, in actual operation the wick must first be 
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primed by lowering the heat pipe to a more nearly level position prior to 
its operation. The advantage of not requiring self-priming is that higher 
maximum heat-transfer rates can be achieved. 

The boundary conditions for the differential equations are 


j P gh w if g i o 

1 4a/Dy $ if g = o 


(7) 


where D 


and 


vs is 
ift(o) 

the hydraulic diameter of the vapor space, 
= M(o) = o , 

(8) 

4>(o) 

- J 

(9) 

P(L) 

= p f . 

(10) 


If the final porosity <j>(L) = is specified rather than P f , then 
is calculated from Eq. (1) with H = 1.955 . 

The differential equations (4) - (6) are integrated repeatedly with 
the Runge-Kutta routine, Each time, the total heat- transfer rate is adjust- 
ed until Eq. (10) is satisfied, or, for a pre-specified porosity distribu- 
tions, until the stress just matches the wick's critical stress anywhere 
in the interval 0 - x - L . 

INPUT 

The input is in Fortran NAMELIST form. The required parameters are 
defined and discussed below. An input form is given in Table 1. 


Headings 

After writing on the first card or line the NAMELIST identifier 
$GRDATA, the user then inputs two lines of descriptive information by 
writing on one card or line HD1 = 60H followed by up to 60 characters 
of title and on the next card or line HD2 = 60H followed by another 
60 characters. GRADE will print these two lines at the beginning of 
the output. 
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TABLE 1 


HD1 - 60H 
HO 2 = 60H 


LIQ = 
TKELVN * 
RHOL = 
ST = 
VIS L = 
HFG = 



GEE = 


INPUT ONLY IF LIQ =0.0 


AW = 

HW = 

diaf = 

PHIO = ) 

PH I F = V HOMOGENEOUS WICK? SET PHI F = PF = 0.0 

pc - \ AND PHIO - WICK POROSITY. (IF ONE SETS PHIO * 0.0 OR 

7 PHIO = DIAF = 0.0, GRADE PROVIDES OPTIMUM VALUES FOR 

DIAVS = THESE.) 


QD0T = 
NQ = 
XQ = 
FQ - 


NELE = 
XELEV = 
ELEV = 


DX = 


NCASE = 

LASTPHT = 

ISLFRRM = 

o Vt- 9 

SEND 
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Fluid Properties 


GRADE automatically computes the required 

fluid properties for one of 

several fluids, which 

the user specifies by selecting a value of LIQ from 

the following list: 




Fluid 

Temperature Range 

LIQ = 1 

Water 

(32F < T < 4 OOF) 

LIQ = 2 

Ammonia 

(-108F < T < 190F) 

LIQ = 3 

Methyl Alcohol 

(-140F < T < 380F) 

LIQ = 4 

FREON-21 

(-55F < T < 305 F) 

LIQ = 5 

Ethane 

(-135F < T < 80F) 

LIQ = 6 

Methane 

(-280F < T < 120F) 

LIQ = 7 - 

Nitrogen 

(-340F < T < 250F) 

The properties are for a temperature TKELVN that the user inputs in 

degrees Kelvin. All « 

other fluid properties are 

automatically computed for 

that temperature. If another fluid is used, set LIQ = 0 . Then, values 

must be specified for 

the following quantities: 


Quantity 

Symbol 

Units 

Liquid density 

RHOL 

Kg/cu. m 

Surface tension 

ST 

Newtons/m 

Liquid viscosity 

VISL 

Newton-sec/sq. m 

Latent heat 

HFG 

Joules/Kg 


Caution must be used with fluids other than ammonia at ambient 
temperatures. Since GRADE was specifically written for ammonia heat 
pipes, the vapor-flow pressure drop was neglected. For fluids with a 
relatively low vapor density, GRADE must be modified to include vapor- 
flow pressure drop. 
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Geometrical parameters 


The following quantities specify the 

cross-sectional geometry 

Quanti ty 

Symbol 

Uni ts 

Wick area 

AW 

sq. cm 

Wick height from tube bottom 

HW 

cm 

Fiber diameter 

DIAF 

cm 

Initial porosity 

PHIO 

dimensionless 

Final porosity 

PHIF 

dimensionless 

Vapor-space hydraulic diameter 

DIAVS 

cm 

Final stress 

PF 

cm liquid 


For non-zero-gravity calculations, if the user sets PHIO = 0.0 , 
GRADE calculates a value just low enough for the wick to self prime at 
the condenser end. If the user chooses to specify the final stress, 

instead of a final porosity, he sets PHIF - 0.0 and PF to a non-zero 
value. Otherwise, PF is set to zero. 

To calculate the performance of a homogeneous wick, both PHIF and 
PF are set to zero, and PHIO is set equal to the desired constant 
porosity of the wick. If PHIO is set to zero, GRADE calculates an optimum 
value that gives the highest capacity for the operating condition. If, in 
addition, the fiber diameter DIAF is set to zero, GRADE will supply 
optimum values for both DIAF and PHIO . 

Heat Input 

The user specifies the heat-input distribution by specifying the 
fraction of the total heat-transfer rate for up to ten segments of the 
heat pipe. Heat transfer is assumed to be constant along the given 
segment. Values for the following parameters are required: 

NQ - the number of segments, which must not exceed ten, into 
which the heat pipe is divided. 
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XQ(I) - the length of the I^ 1 segment in cm. The segments must 
be numbered consecutively along the heat pipe beginning 
at the condenser end. 

h 

FQ(I) - the fraction of the total heat-transfer rate entering the I 
segment. If the segment is adiabatic, FQ( I ) is 

negative. 

QDtf)T - A nominal heat- transfer rate in watts, which is the user's 
best guess at the maximum. A close guess reduces the 
number of iterations to the final answer. 

Elevations 

The user specifies the heat-pipe orientation in a gravitational 
field by inputing values for elevations of points along the heat pipe 
where the slope changes. Between points, GRADE assumes a linear variation 
of elevation. Values of the following parameters are required (except 

for zero gravity): 

NELEV - the total number of points along the heat pipe, which must 
not exceed 10. 

t h 

XELEV(I) - the distance along the heat pipe in cm to the I point. 

Both ends of the heat pipe must be input, therefore the first 
point must be at zero distance [XELEV(l) = 0.0], and the last 
at the total heat-pipe length [ XELEV( NELEV) = L]. 

ELEV(I) - the elevation in cm of the 1^ point relative to a horizon- 
tal refernece plane. 

GEE - gravitational acceleration in standard gravities. 


Program Control 


The following quantities relate to the mechanics of the program and 
user options: 
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DX - integration step size, in cm. 

ISLFPRM - equals 1 if the user requires the wick to self-prime at the 
specified elevations and zero load. 

- equals 0 if the user does not require the wick to self prime. 

NCASE - equals 1 if another case follows. 

- equals 0 if another case does not follow. 

LASTPHI - equals 1 if the porosity distribution of the last case is 

to be used. 

- equals 0 if it is not to be used. 

$END - ends NAMELIST input. 

If another case is to follow (NCASE =1), then another NAMELIST input 
follows that is exactly like the first except only those quantities that 
are to be different in the new case are included. Up to ten cases can 
be run. The NAMELIST input for the sample calculations are included in 
Appendix A. 

OUTPUT 

GRADE first prints the headings that the user has input. The input 
variables and computed thermodynamic properties are printed next. The 
final iterated solution follows, which consists of the maximum heat- 
transfer rate, the total amount of liquid in the wick, and a tabulation 
of distance along the heat pipe, stress, static head, porosity and critical 
stress. The static head is that part of the stress due to gravity. The 
critical stress is the maximum the wick can sustain without failing. 

SAMPLE PROBLEM 

To demonstrate GRADE, a typical heat pipe was analyzed. It has a 
,356-cm x 1.128-cm (.140-in x .444-in) slab wick with a fiber diameter of 
.01016 cm (.004 in). The heat pipe is 180-cm long with single condenser. 
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adiabatic and evaporator sections of 60 cm. Performance was studied for 
ammonia at 294°K (70° F). 

For the first run (output in Appendix A) the evaporator end was 
elevated 2 cm. A final stress of 7 cm is specified, which is slightly 
above the pumping ability of a ,0178-cm (.007-in) pore. (We have previous- 
ly judged that such a pore corresponds to the critical stress at which 
grooves lose communication with the wick. By setting PHI 0=0. we let 
GRADE calculate the initial porosity that is just low enough to prime the 
condenser end of the wick. We require in the first run that the entire 
wick has a no-load self-priming capability at the 2-cm elevation (ISLFPRM - 
1). GRADE computed a maximum heat-transfer rate of 81.3 watts (3843 watt- 
in). The stress and porosity distrubtion are shown in Fig. 1. Until the 
break in the stress curve at 120 cm, the reduction in porosity is due to 
the self-priming requirement. 

For the second run, the same porosity variation was studied in zero 
gravity. The calculated maximum heat-transfer rate is 141 watts (6661 
watts-i'n). Notice that in zero gravity, the initial stress has the low 
value of the vapor-space capillary pressure, whereas in earth gravity the 
initial stress is set by the hydrostatic head at the top of the wick 
of 1 .128 cm. 

For the third run, the self-priming requirement of the first run 
was dropped (ISLFPRM = 0). This increased the maximum heat-transfer rate 
to 115 watts (5452 watt-in). The wick porosity remains constant at its 
initial value until at 62 cm the stress reaches the capi 11 ary-pressure 
limit of the wick for emptying. 

In the fourth run, the wick designed in the third run is studied 
in zero gravity. The maximum heat transfer rate is 200 watts (9439 watt-in) 

In the fifth run, we study a homogeneous wick with a porosity of 
.683, which is the same final value as the graded-porosi ty wicks. Thus 
the homogeneous and graded-porosi ty wicks have the same elevation at which 
the heat- transfer capacity goes to zero. The calculated capacity for the 
homogeneous wick at a 2-cm elevation is 30 watts (1409 watt-in). 
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STRESS (cm liquid) 



DISTANCE (cm) 

FIGURE 1. Porosity and Stress Distribution for the Sample Problem 
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APPENDIX A 


SAMPLE PROBLEM INPUT AND OUTPUT 
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NAfCLIST INPUT 


iGROATA 

HO 1=6 0 HORACE SAMPLE PRO IL EM 
HO2=6QHSELF-PRIHIN0 WICK 2-QM ELEVATION 

ua=2 

TKELVN=294. 

GEE-1 • 0 
AW = 0* 401 
HW =1 • 123 
DIAF=a,J101E 
RHIC=0.0 
PH IF=0 • 3 
PF =7. 3 
DI AVS = G • 47L 
QOOI=10Q. 

NO = 3 

F Q -*• i.C f C.u»i,0 
XQ= oj « *60* ■» 6 G • 

NELEV = 2 

XE LEV = 0 • 0 * 18 H • 

EL £V= 0 • ‘J « 2* C 
DX =2* 

LASTP HI = Q 
IS LFPRM=i 
NO ASE = 1 
SEND 
fGROATA 

HO 2=6 OH SSL*-' -PRIMING WICK ZERO GRAVITY 

GEE=C . G 

LA$rPHI=l 

SEND 

SGRDATA 

HO 2=6 OH iVON-StL^ -PRIMING WICK 2-CM ELEVATION 

GEE = 1 « 3 

IASTPHI=3 

IS L r F RM = J 

SEND 

SG R3A T A 

HO 2=6 3HN3N- St L^-P RIMING WICK ZERO GRAVITY 

GEE=u . 0 

LA STPHT = 1 

if NO 

SGRDATA 



ri02^blHr»OM0Gf.M£aJS HICK 2-Ctt ELEVATION 

GtE-1 • 0 

LA STPHI =Q 

PF=0, 3 

PHlQ=.6iT5 

NCASE=C 

I3LFPRH=3 

ftE NO 



GRADE SAMPL C PROBLEM 
SELF-PRIMING 41 CK 2-CM ELEVATION 

INPUT VARIABLES ANO ^LJIB PROPERTIES 


LIQUID NUMBER 

TEMPERATURE, . . . 

LIQUIJ DENSITY ..... 

SURFACE TENSION. 

LIQUID VISCOSITY 

LATENT HEAT... 

GRAVITATIONAL ACCELERATION. . . . 

HICK AREA 

HICK HEIGHT . 

HICK c I3ER DIAMETER 

INITIAL POROSITY 

FINAL POROSITY 

FINAL STRESS.. 


VrtPQR- SPACE DIAMETER, 


NOMINAL HEAT-TRANSFER RATE. 
NO. HEAT-INPUT SECTIONS.... 
SECTION NUMBER 1 

SECTION LENGTH 

Hr A T -IN T FRACF ION... 
SECTION NUMBER ?. 

SECTION LENGTH 

H CAT-IN PUT FRACTION... 
SECTION NUMBER 1 

SECTION LENGTH 

HEAT-INP'JT FRACTION... 


NO. ELEVATION POINTS..,. 

ELtVATION POINT NO, 1 

J I ^ T ANC i TJ C I N T 

elevation jp foi n r 


CASE NO. 1 


LIQ = 2 

TKELVN = 2. 943 0 G £«■ G 2 ) EG REES KELVIN 

RHQL = b,G9C35t>D2 KG/CU. M 

ST = 2 . 0 94 0 9 E> A 2 NEWT ONS/M 

VI SL = 1.37440E-G4 Nt WTUN-St C/SQ* M 

HFG .= 1.1A576EKJ6 JOUL c S/KG 


GEE 


1 -. aeacct+oc- 

STANJARO 

AW 

z: 

L.ulDCCE-01 

SO. CM 

HW 

- 

i*iaaoot>aP 

CM 

01 A r 

- 

l.Ui&0Cc-T2 

CM 

D H 1 0 

- 

a. 


PHIF 

— 

p . 


PF 

- 

7. DOaC'CLfOO 

CM LIQ. 

Ut AVS 

= 

4.700CCE-31 

CM 

i Li 0 T 

- 

I * y G u L L E + j 2 

WATTS 

NQ 

~ 

Tt 


XQ 

- 

6.coooge*ci 

CM 

FQ 

— 

-i,0'jaC'OE+cj 


XQ 

= 

t * Q 0 j 0 t. E + C 1 

0 '1 

FQ 

= 

c. 


XQ 

— 

6 . 0 r } G 0 R + -J 1 

CM 

FO 

z. 

1 * CioO-jOL + yP 


NELCV 

= 

o 

(- 


<~Lr;.V 

=r 

r* 

« 

' M 

r t. E V 

= 

p 

!-■ • 

C M 



ELEVATION POINT NO. 2 

OTST ANC£ TO POINT X€LE V = 

ELEVATION OF POINT HLFV = 

INTEGRATION STfcP SIZE DX = 

ANOTHER CASE (C = NO, i=YESI.*.. NC ASF = 

USE LAST POROSITY OISTN LA5TPHI = 

REQUIRE SELF PRIME ....... * ... . ISLFFRM = 


1 . a G 0 n C E ♦ G 2 EM- 
,2.0a0GL:+aG CM 


2. G03 OOE«-UO CM 

1 

r 

v- 
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FINAL SOLUTION 


THE MAXIMUM HEAT-TRANSFER RATE 


TOTAL LIQUIO 

IN WICK IS 

OISTANCE 

STRESS 

(CM ) 

(CM LIQ.) 


0. 

l.iesaE+oo 

2*0aGGE+0D 

1.1505E+00 

4.0000E+0C 

1.1737E+00 

6.QJQ0C+UG 

1 • 19 75 E+ Q 0 

0.OQOQC+OO 

1*2220 E+-GQ 

l.OGOOF+Ci 

1.2473E+0Q 

1 • 2 0 0 OE + G 1 

1 • 27 34E + 0 0 

i*40.aGE + Gl 

1 • 300 4E + G 0 

1 .eooaE+oi 

1 • 3282E+Q 3 

i.aaouE+oi 

1.3570E+30 

2* 09G0E+C 1 

1 .3360E+TO 

2.20C0E+G1 

1.4176E+00 

2 • 4QG G£ + Q 1 

1 • 4496E+03 

2.60G0E+C1 

1.4327E+0C 

2.0QOCF + O.1 

1 • 51 7 G E+ 30 

I. OGUOE+Ci 

1 • 5525E+ OQ 

J.20QQE+ 01 

1. 5994E>aa 

3.43QuE+Cl 

1* 6277E+0G 

3.60G0E+L1 

1.6674E+DQ 

3.80G0E+G1 

1 • 7 u 06 E+ Q Q 

4. OOCGE+Ql 

1. 7514E+G0 

4*200 QE + u 1 

1.7 *532*01 

4*4QU0t+9i 

1*841 9E +03 

4 • 6 QO 0 1 + G 1 

1.3398t+u0 

4.890'QC + Ci 

i* 9 S95E + 03 

5*Q0COt.+Cl 

1. 9912E+ 03 

5.2JG0E+01 

2* 0 44 3 E + G 0 

5 . 4 0 0 Q E + 0 1 

2* 1004E+UO 

5.6000E+01 

2*1 582 £+00 

5. 8000 E +01 

2 * 211 lfc>oa 


3 8* 14453E+01 

WATTS 


433Q6E+01 GRAMS 



STATIC HEAD 

P0R03I TV 

CRIT* STRESS 

(CM LIQ*f 


(CM LIQ.) 

1*1200 £ + 0 0 

8. 7265E-C1 

2*2 052E + 0 0 

1.1E32E+00 

0. 7046E-31 

2.2487E+00 

1.1724E+GC 

8. 6829E- 01 . 

2.2 921E+G0 

1* 1947E+QQ 

6.6613E-U 

2*335 6E +00 

1.2169E+0G 

0* 639 3E-G 1 

2.3790E+90 

1.2391E+GC 

8.6104E-G1 

2* 4225E+0 0 

1 « 2613E + CG 

8.5971E-G1 

2.4659E+00 

1.2836E+0G 

8* 5759E- G 1 

2*509 4E +00 

1.3058E+00 

8* 5548E-01 

2 • 5 52 8E + 0 0 

1.3230E+GQ 

8* 53 3 QE- 01 

2* 5962E+0 C 

1 ♦ 35C2E+ 0 Q 

3. 5129E-01 

2.6397E+00 

1.3724E+00 

8. 492 lL-01 

2.6831E+0C 

1* 3947E+Q0 

0. 4714E-Q1 

2.7266E+00 

1.4169E+0D 

3 * .4 5 J 9 £ - u i 

2*770 OE + 30 

1. 4391E+00 

8*43 u 4E- 01 

2*013 5E+0C 

1 • 4613E+0C 

8*411 GE-G1 

2*856 9t +00 

1.4836F +00 

8. 3097E-Q1 

2* 90 0 4E + 0 G 

1*50 5dE + QU 

8. 3695E-01 

2.9430E+Q 0 

1*520 GE+OG 

8. 3494E-01 

2*98 72 £4-00 

1.5502E+0Q 

8* 32 9**E-0 1 

3. T307E+3 0 

1.5724E+0G 

8 • 3 u 95E-0 1 ' 

3 • 0 74 1 E + 0 0 

1. 5947E+0U 

8. 2897E-01 

3*117 6E + C C 

1.&169E+C0 

8. 2 70 Ot * 01 

3*161 OE+OC 

1*6 391 1 ♦ 0 0 

8. 25C4C-01 

3*2 045E+ 0 0 

1.6613E+0C 

8. 23u9E-Gi 

3.2479E+QD 

1.6836E + 00 

ft. 2U4E-31 

3. 2914E+00 

1* 7O50E + GG 

(». 1921E-01 

3*334 0f +00 

1*7200 6. + CO 

ft. 1729E-G1 

3.3732E+0C 

1.7532F+00 

ft. 1537E-P1 

3 *421 7L +00 

i. 7724E+0C 

8.1 346E-01 

3.4651E+0C 





6.03005+01 

2.28945+00 

6.20005+01 

2. 3441E+QG 

6.40005+01 

2 .4087E + UQ 

6.60025+01 

2. 47445+0 0 

6.80035+01 

2.5411E+00 

7.03005*01 

2 . 6u 88E + 0 3 

7.2 JuUE+Ci 

2* 67 765 + 9 u 

7 ,40036+01 

2. 74745+03 

7 *61005 + 01 

2. 61335+03 

7.90035+01 

2. 39035+00 

8. COQOE + Oi 

2. 96345+00 

3* 2000 £ + 01 

3*03765+00 

3.4ao0£+01 

3* 11305+03 

3*6009.5 + 01 

3*13455+00 

3. 33095 + 01 

3 « 25 7 25+ 0 0 

9.00035+01 

3 • 34615 + 0 0 

9.23005+01 

3*42625+00 

9,43005+01 

3.50755+00 

9.63CCE+01 

3 » 590 0 E+ 0 0 

9,30005+01 

3.67335+30 

1. GO 005 + 02 

3. 75395+30 

1.G2GGE+.0 2 

3 . 84525+00 

1 , 04005+02 

3 * 932 35 + 0 0 

1*06005+02 

4.921 3 E + uQ 

1. 03005 + 02 

4. 1121E + 0 0 

1,19005+02 

+ . 2Q37E + 00 

1*12005+02 

+ * 2 96 75 + 00 

1.14095+02 

4. 39115+30 

1 *16005+02 

4.43685+00 

1.U33E + 02 

*4,53 435 + 03 

1.20 9 0F + C2 

4.63275+30 

1.22005+02 

4.73196+00 

1.24005+02 

4.87995+ j3 

1,26095+02 

4.97706+30 

l,23Gu5+C 2 

5. 07425+60 

1,33005+02 

5. 171 75+ JO 

i . 32005 + 02 

5, 2 6916 + u 9 

1.3+035+02 

5 • 3 6 6 4 E + 0 0 

1 • 3 5 0 0 1 + 02 

5 .46335+30 

1*33005+02 

5.55975+30 

1.4JG05+C2 

5 . 65535+23 

1.42005+02 

5 . 75G 05+ CO 


1.7 947 c + C 0 

8.11 565- G 1 

1. 3169E+0Q 

0. 0 9635-3.1 

1*83 91 E+ 00 

8. 0 7836-01 

1.86135+00 

8, 0592 £-01 

1*88365+00 

8. 0406E-G1 

1.9G58E+GQ 

8.U221E-01 

1.9250E+00 

3. Q0 365-01 

1*95025+00 

7. 98525-01 

1*97245+00 

7. 96695-01 

1.9947F+CC 

7*94375-01 

2.0169E+Q3 

7.93066-31 

2. 03915+00 

7. 91265 -01 

2. 0613E+CD 

7. 89466-01 

2* Q836E + CG 

7. 87675-01 

2.1156E+C0 

7. 85895-01 

2.128UF+0G 

7. 84125-01 

2 • 15025 + a 0 

7. 32 366-01 

2* 17245+00 

7. 80635-01 

2* 19475 + 00 

7. 78355-01 

2.21695+00 

7. 77115-Ci 

2* 2 3915 + 00 

7. 753 35- u 1 

2» 2613 5 + 00 

7. 73665- Cl 

2.28 36 E + 00‘ 

7. 7194E-C1 

2.305 35+00 

7,70235-01 

2. 323CE+C0 

7*68535- Cl 

2,35026+00 

7. 66835-01 

2. 57245+00 

7, 6515 E - 0 1 

2. 39475+00 

7, 63476-01 

2 ■ 4169F+0G 

7.61805-01 

2. 43915+00 

7.60135-01 

2,46135+00 

7. 53475-ul 

2. 48365+00 

7. 5 632 z. - Cl 

2. 50536+0 3 

7. 551 3c- 01 

2.52805+03 

7 * 5224 5 — Cl 

2. 55026 + L 0 

7.48625-01 

2457245+00 

7. 45025-01 

2.5947E+0Q 

7*414oc-C1 

2, 61695+00 

7. 379+5-01 

2.63 91 E + 0 0 

7 • 34465-01 

2.66135+00 

7. 51 J 3E-G1 

2.68365+00 

7, 27675- Cl 

2 * 7 0 5 3 5 + 0 u 

7, 2 4 3 6 5 - 0 1 


3.508 65 +0P 
3. 552 Ofc + 0 C 
3.5955E+uC 
3* 6389E+0 0 
3*6 82 45 + 0 0 
3.725 85 +3C 
3*76925+00 
3,31275 +3 C 
3*856.15 + 00 
3*89965+00 
3*94305+00 
3*93655+0 C 
4*0 2995 + 00 
4,37345+OC 
4.11685+00 
4.16025+ JC 
4*20375+00 
4.24715+00 
■4*29 365+0 0 

4 * 3 3 4 G E + 0 0 
4,3775E+G C 
4. 42096+0 C 
4.4644E+C0 
4. 5 L 785+ 3 0 
4,55126+0 C 
4*59475+00 
4*63816+0 G 
4.68165+30 
4*72506+00 
4*76855+00 
4.31195+00 
4* 3 5 5 4 5 + G C 
4.8)946+30 
4*97705+0 0 
5* C 742“ + jC 
5. 1 7175+0 C 
5*26915+11 
5. 3 6 o 4 - + C 0 

5 • 4 6 3 ,’> r + f. ; 0 
5.5597 r. ♦ j C 

5 • 6 5 5 3 6 + 0 C 

6 • 7 5 0 f E + 0 0 



1.44Q 36 + 0 2 
1.4S&QE+02 
1 • 45 OQE + 0 2 
1.53QQE+Q2 
1.5200E+02 
1.5403E+C2 
1.5600E+G2 
1 • 58 0 0E + 0 2 
1.&0D0E+02 
1 ♦ 6 2 o 0E + C 2 
1 * 64 0 0E + O 2 
1.6600E+02 
1.6300E+02 
1.7300E+02 
1.72G0E+02 
1.74096+02 
1. 7600E+02 
1.73006+02 
1.8300E+G2 


ro 

o 


5.8435E+GQ 

2*7290 E + C 0 

5 *93556+03 

2. 7532E+C0 

6.02596+30 

2 • 7 72 46 +00 

6 .11416+00 

2.7947E+00 

6.20016+09 

2.3 169 E + 0 0 

6.23366+00 

2.83916+00 

o • 36426+ 0 0 

2.9613E+00 

6.44166+00 

2.8836E+C0 

6.5157E+00 

2.9053E+00 

6.5360 6 + 0 0 

2.9290E+CG 

b • 6522E + u 0 

2.95Q2E+CG 

&. 7141E+00 

2.97246+00 

6.7 7 156+30 

2.99476+C0 

6.32 406+ 0 0 

3. 01 69 £+00 

6.37146+00 

3.0 391 E + G 0 

6.913 56 +00 

3.3 6136+0 C 

6. 95016+00 

3 * 0 636 E + 0 D 

6 • 9 31 1E+ 0 0 

3. 10536+00 

7 • 0 0 62t+ 00 

3.129GE+0U 


7. 2113E-01 

5 • 3435E + 0 0 

7.179 3E-0.1 

5.935 5E +00 

7 • 14916-01 

6. 02586+00 

7. 1194E-01 

6.1141E+0G 

7. Q997E-G1 

6 • 2 GO1E + 0 G 

7.06306-01 

6.2836E+0C 

7. U365E-U1 

6.3642E+UC 

7.011 2E— G 1 

6.441 6£ +00 

6.9872 E -01 

6 • 515 7£ + 0 0 

6. 9646i£-C 1 

6.5860E+OG 

6. 94346-01 

6.6522E+0P 

6. 92 3 7E-C 1 

6.71416+00 

6. 9C55E-01 

6.771 5t + 9 0 

6. 889GE-C1 

6.82406+00 

6. 3741E-01 

6.8714E+0C 

6. 96106-01 

6. 91356 + 00 

6. 64 9 66- Cl 

6.9501E+00 

6.34006-31 

6.99116*00 

6. 8 322 E- Cl 

7.00626+30 


i 



GRADE SAMPLE PROBLEM 
SELF-PRIMING wick zero GRAVITY 


INPUT VARIABLES AND FLU 10 PROPERTIES! 


LIQUID NUMBER* * * * 

TEMPERATURE 

LIQUID DENSITY. 

.SURFACE TENSION.. 

LIQUID VISCOSITY 

LATENT HEAT................... 

GRAVITATIONAL ACCELERATION.... 

WICK AREA. 

WICK HEIGHT 

WICK FIBER DIAMETER * 

INITIAL POROSITY 

FINAL POROSITY... 

FINAL STRESS 

VAPOR-SPACE DIAMETER 

NOMINAL HEAT-TRANSFER RATE.... 

NO. HEAT-INPUT SECTIONS 

SECTION NUMBER 1 

SECTION LENGTH 

HEAT-INPUT FRACTION...... 

SECTION NUMBER Z 

SECTION LENGTH 

HEAT-INPUT FRACTION...... 

SECTION NUMBER 1 

SECTION LENGTH.. 

HEAT-INPUT FRACTION 

NO. ELEVATION POINTS.... 

elevation point no. i 

QISTANCE TO POINT 

ELEVATION OF POINT 


CASE MO. 2 


LIO 

— 

2 


TKELVN 

' = 

2*940 CGF*G2 

DEGREES KELVIN 

RHQL 

= 

6. Q9435F+ 02 

KG/CU. M 

ST 

= 

2 • 0 94 G9£- 0 2 

NEWTON3/M 

VISL 

r 

1 « 37840E- 04 

NEWTON-SEC/SQ. 

HFG . 

- 

1.18576LH06 

JOULES/KG 

GEE 

= 

0 . 

STANDARD GRAVI 

AW 

- 

L.OIOGOE-Dl 

SQ. CM 

HW 


1.128G0c+00 

CM 

DIAF 

r 

1.G164GE-Q2 

CM 

PHIG 


G. 


PH IF 


C. 


PF 

= 

7.00GSDE4-0G 

CM LIQ. 

DIAVS 

= 

4* 70 0 G CE- 0 1 

CM 

COOT 

- 

1.G00UCE+C2 

WATTS 

NQ 

~ 

3 


XO 

s 

6. D000CE> 31 

CM 

fq 

= 

-l.COaCoE+G'; 


XO 

= 

6 • G 0 0 G CE* 0 1 

CM 

FQ 

- 

0. 

- 

XQ 


6. 0JGGCF + 31 

CM 

FQ 

= 

l.GOQGOE+jn 


NELEV 

= 

2 


X£ LEV 

r 

0. 

CM 

ELEV 

5 

G. 

C* 



ELEVATION! POINT NO. 2 

DISTANCE TO POINT . XELEV = 

ELEVATION OF POINT . EL£V = 

INTEGRATION STEP SIZE.......** OX = 


1.&D0D0E+C2 
2. C jOOOE^OO 

2.0QODCE4-OO 


ANOTHER CASE <0=NO, i=YE$>.... 

USE LAST POROSITY DISTN 

REQUIRE SELF PRIME 


NCASE ■= 1 

LASTPHT = i 
ISL C FRM = 1 



final solution 


THE MAXIMUM HEAT-TRANSFER RATE IS 1.39551E+02 WATTS 
THE TOTAL LIQUID IN WICK IS 3.433Q6E+01 GRAMS 


DISTANCE 

(CM) 


0. 

2*Q900E + iiG 
4. JUJOE+GG 
6 * 0 J u Gt +Q 0 
3* G 0 0 9 E + i> G 
i.0aG3E+01 
l * 2 a o je + o 1 
1.4009E+G1 
^ 1*6900 E +01 

w l.aao DE+oi 

2*0000 £ + G 1 
2.20C0E+01 
2.40Q0E+C1 
2 # 69 0 Ct + 0 1 
2.6OOOE+01 
3*0090 E +01 
3.23001+01 
3.43<JuE + Gl 
3,6aOQf+Gi 
3*60001+01 
4* Q0O3E + C1 
4.23C0E+U1 
4*43031+01 
4*60001+01 
4.3QL0E+C1 
5.09CCE+01 
5.23OOE+01 
5*40001+01 
5.60C0E+O1 
5*30036+01 


STRESS 

STATIC HEAD 

(CM LIQ* > 

(CM LIQ.) 

2 .9840E-01 

0. 

2 • 93921* 0 1 

0* 

3*095 it- 3 1 

a* 

3 • 0325E-01 

9. 

3.0722E-01 

9. 

3*124 sHE- 01 

0. 

3* 1912 E- 01 

a . 

3*272 IE- 0 1 

0. 

3*36851-01 

0. 

3*4310 E- 01 

0* 

3*610 6E- 0 1 

a. 

3. 7582E- 0 1 

a. 

3 • 9246E-0 1 

n. 

4* 11G SE-G 1 

0. 

4 • 31 7 3 E- 01 

0. 

4* 5464E-01 

0. 

4 • 7977E- 0 1 

0* 

5* 0727 E- 01 

3* 

5 • 3 724E-G 1 

3. 

5. 63791-01 

0 • 

b. C503E-G1 

c* 

6*4 3 0 71-01 

0. 

6*840 IE- 0 1 

a. 

7*2 79SE-01 

o • 

7.750 9E-bi 

0. 

3 • 254 7E-0 i 

0* 

3. 7922E-01 

3* 

9.3649E-91 

Q. 

9*974 0E- 0 1 

3 . 

1.0621E+Q9 

0 * 


PQROSITf 

CRIT. STRESS 


(CM LIQ* 1 

fl. 7265E-91 

2.2C52E+00 

8. 7046E-&1 

2*24871+00 

3 • 68291-01 

2.2921E+J0 

8.6613E-C1 

2* 3 3 56E + 0 P 

8.6398E-01 

2.37901+90 

8. 61341-01 

2*422 5E+0G 

8*597 IE- 0 1 

2*465 9£ + 0 0 

8, 5759E-01 

2.5J94E+CO 

8* 5548E-Q1 

2* 5 52 8F+ 9 0 

3.5333E-91 

2 • 59621+0 n 

8* 51231-CI 

2*63971 + 00 

8.492IE-U1 

2.6331E+9C 

6*471 41-01 

2*72 &6E + 00 

6. 45091-01 

2*770 Cl + 00 

8. 4 3 04 E -Cl 

2*813 5t + 9 P 

e.4iooE-9i 

2 . 8 5691 + 0 n 

8. 3397E-C1 

2.9004^ + 00 

8. 3695 E- 01 

2*94381+00 

8 • 3 4 9 4 E - 0 1 

2.93721+Cn 

6. 3294E-P1 

3*03071 +00 

8. 3095E-H1 • 

3*0741 E + GO 

8. 23971-61 

3.11761 +00 

8* 2 7b OE -u 1 

3. 16101 + 0 0 

8. 2504L-L1 

3*2^451 + 0 G 

B. 23G9E-C3 

3*247 9P +00 

3. 2U4E-D1 

3. 29141+0 C 

3* 19212-01 

3* 334 3E + 0 0 

8. 1729 E -01 

' 3. 37821+00 

8. 1537E-G1 

3*421 75 +00 

8*1 J4&1-G 1 

3*46511+00 



6.0:]QI3£ + ul 1*130 66 + G 0 
6.23006+01 1 • 23 1 7E+ G 3 

0. 40006+Gi 1.27446+00 
6.63GQE+01 1.3483E+G0 
6.30006+01 1 • 425 3E + G 0 
7.30006+01 1.53296+33 
7.2 GGQ6+Q1 1 • 5527E+uO 
7,4CC0£+ul 1 • 664 3E+ GO 
7.6 JJOE + Cl 1 • 74 772+ C 0 
7 • 8 0 J C £ + 0 1 l*33332+yQ 
3.0Q0OE+G1 1*92026+03 
8.2003E+G1 2*0093 £+0 3 
3 *4 J 006 + 0 1 2* 100 4E + 0 0 
3.60032+01 2. 19346+00 
8.80C32+01 2*2385 E+ 00 
9.QU00E+C1 2*33552+00 
9.20002+01 2.4347E+00 
9*40 00 E + 01 2*53592 + 00 
9.63CGE+01 2*63932+03 
9 . 3 O 0 OF+QI 2. 79476+30 
1 • 3 0 0 0 £ + 0 2 2.9324E+00 

K 1. 02002+02 3. 01222+00 
1*04302+02 3*12436+00 

1. G6CGE+G2 3*23366+00 
1.03006+02 3.3552E+00 
1 . 1QUQE+0? 3.47422+00 
1.1200E+0? 3*59546+ QQ 
1.1400 £+32 3. 71912+03 
1.1603E+G2 3. 34516+00 
1.18002 + 32 3* 9735E+03 
1.23002+02 4.10456+00 
1.22CCE+02 4.2364E+G0 
1.2 4002 + 02 4. 36626 + G 3 
1.26036+02 4.49436+30 
1.2300E+C2 4*62346+00 
1*33002+02 4. 75236+00 
1* 320 06 + G 2 4 • 6312E + 0Q 
1.3+GJE+02 5. CO 3 72 + 03 
i. 3603 £+02 6*13776+00 
1. 330 02+0 2 5. 23482+03 
1. 43 C 06+02 5.39366+00 
1.42006+02 5.51472+03 


o a o q a a -j o c o u a o a o j o a c o o o o □ a 


fl 

V • 

8.11562-01 

3*50366+00 

3. 

3.G963E-G1 

3.55206+ nr 

3. 

8. C7302-G1 

3.59552+00 

0 . 

8*35922-01 

3 . 6 3 8 9F + u C 

a. 

8.34062-01 

3.68246+00 

0. 

8*22212- 01 

3*725 8 E+0C 

0. 

8. 00 366- Cl 

3*76926+00 

j . 

7, 9352E-Q1 

3.312 72 +0G 

0. 

7, 966 9E-ij1 

3*85616+00 

0. 

7. 94876-01 

3.8996E + C C 

3. 

7. 930 62- Cl 

3.94306+00 

0. 

7. 91262-31 

3.98656+00 

0 . 

7. 39462-01 

4.02 996 + C 0 

3. 

7. 3 7 6 7 2 - G 1 

4. Q 7346+0 0 

0. 

7, 35896-01 

4.1168E + QC! 

3 . 

7. 64122-01 

4 . 1602E + 0 Q 

0. 

7 . 3236E-Q1 

4.2G37E+0C 

a. 

7* 8G606-C1 

4* 24716+ 0 C 

o . 

7. 78 252-0 1 

4.29062+00 

0 • 

7.77112-01 

4.334CE+0C 

0. 

7 . 7 5782-ul 

4. 37756+0 G 

0 . 

7. 73662-01 

4.42092 + 0 0 

0. 

7.71942-31 

4.46446 + 0 C. 

o. 

7 . 7323 ti-Gl 

4*53736+00 

0. 

7. 68532-01 

4.5512F+0C 

Q. 

7. 5633I-Q1 

4.59472+ n 0 

a. 

7 . 6515 2 -Cl 

4,63812+00 

0. 

7. 6347E-ul 

4*68162+0 G 

3. 

7* 613 32-01 

4.77506 + 0.0 

0. 

7.6013E-C1 

4. 76352 + 0 G 

a. 

7 . 58472-01 

4*81196+0 C 

0 . 

7. 56622-01 

4.3 5546 + 0 0 

3 • 

7. 55132-C1 

4 • 8 9 8 8 F +Q0 

0. 

7*52242-01 

4.977G2+0C 

1. 

7.48622-01 

5. a 7422 + 9 0 

G . 

7. 450 2 2 -Cl 

5.1717E+jC 

G. 

7.41462-Q1 

5.26916+f ? 

3. 

7, 37942-ri 

5 * 36646 + 0 C 

0. 

7. 3 446 E-.G 1 

5.46336+0 C 

3 . 

7. 51 w 3 2- Cl 

5 . 55972 + 0 0 

n 

j « 

7. 27672-01 

5. 65576 + 0 0 

3 . 

7 . 24362- 0 1 

6.75006+00 


1 .44 00 E 4-0 2 5.6363E4-00 

1.450GE4-02 5.7564E4-Q3 

1.430QE+02 5.873UE*00 

1.530064-02 5 .9362£tQ3 

1.52GO £4-02 6.0955 E+Gtf 

1.54C0E+C2 6.2QJ)4E+0a 

l#5oDiJF + 02 6.33Q5E4-0O 

' 1.53GQE+Q2 6 . 3 9 5 1 E 4- Q Q 

1 #6u C OE + O 2 6.433 9E 4-0 3 

1.62G3E+G2 6. 5652E4-QQ 

1.64GQE4-02 6.64Ibc4‘€0 

1*65 ODE +C2 6. 7097E+Q0 

1.63Q0E4-D2 b. 7693E* GO 

i.raaoE+02 6.82i7£+ao 

1.72G3E+G2 6.8649E+00 

1.740 3E4-02 6.B939E+GQ 

1.76C0E+02 6. 923 6E 4- 3 0 

1.7300E+C2 6«9385£*3D 

1.3du3E+C2 6 * 9436E4-2U 


ro 

Cn 


ca a o ta o 


a . 

7 » 2113E-01 

5.8435E+3C 

0. 

7 . 179BE-01 

5.9355E+0G 

0. 

7. 1491E-01 

6.025BE+QC 

0. 

7. 1194E-ul 

6* 1 1 41 E 4-0 0 

0. 

7. 3907E-01 

6,2GGlt >00 

0. 

7. 0630 £-01 

6.2836E4-C0 

G . 

7.036SE-Q1 

6.36426*00 

0. 

7. 0112E-01 

6.4416E4-0 0 

0. 

6. 9872E-C1 

6, 5 15 7E* 0 0 

Q. 

b . 964 6 £ - 0 1 

6.5B6GE4-0 0 

a. 

6.9434E-01 

6.6522E4-00 

0. 

6. 9237 £-01 

6. 7141E4-0 0 

0. 

6.90 55 E-u 1 

6.7715E+3C 

0. 

6. 3890E-01 

6. 3240E4-0 G 

0. 

6.3741E-G1 

6.871 4E 4-90 

0. 

6. 86 1 0 £-0 1 

6.9135F+0C 

n. 

fc. 6496E-C1 

6. 95 9 IE 4- GO 

0. 

6. 8400F-Q1 

6.9811 £4-30 

0. 

6.6322£-Ql 

7.0 062 c. 4-00 



2- CM ELEVATION 


CASE NO. 3 


GRADE SAMPLE PROBLEM 
NON- SELF-PR I MING NICK 


INPUT WAR I A3LES AND FLU IG PROPERTIES: 


LIQUID NUMBER 

TEMPERATURE 

LIQUIO OENSITY. 

SURFACE TENSION 

LIQUID VISCOSITV 

LATENT HEAT.. 

gravitational acceleration.. . . 

WICK AREA.* 

WICK HEIGHT 

WICK F I3ER DIAMETER . 

INITIAL P ORO 3 1 T V 

FINAL POROSITY ... . 

ro 

Ch 

FINAL STRESS 

VAPOR-SPACE DIAMETER 

\ 

NOMINAL HEAT-TRANSFER RATE.... 

NO. HEAT-INFUT SECTIONS....... 

SECTION NUMBER 1 

SECTION LENGTH..... 

HEAT-INFUT FRACTION...... 

SECTION NUMBER 2 

SECTION LENGTH 

HEAT-INFUT FRACTION 

SECTION NUMBER 3 

SECTION LENGTH 

HEAT-INPUT FRACTION 

NO. ELEVATION POINTS 

ELEVATION POINT NO. 1 

□ ISTANCE TO POINT. 

ELEVATION OF POINT,. 


LIO 

= 

2 


TKELVN 

= 

2.9LQCGE+Q2 

DEGREES KELVIN 

RHQL 

■= 

6.0 3435E+02 

K.G/CU. M 

ST 

— 

2.G94C9E-02 

NEWTONS/M 

VI SL 


1.3784GE-04 

NEWTON-SEC/SQ. M 

HFG 

— 

i. 13576E+ 96 

JOULES/KG 

GEE 

= 

l.DOQOtt + OG 

STANDARD GRAVITIES 

A H 

- 

4. OlOOOE-Qi 

SO. CM 

HW 

- 

1 * 128 Li GE+ G Q 

CM 

DIAF 

- 

1.01600E-02 

CM 

PH I 0 

- 

0. 


PHIF 

~ 

0. 


PF 

- 

7.000QGE+QC 

cm Liq. 

01 A VS 


4.7000QE-01 

CM 

QDQT 

= 

1.0Q0CC E + 02 

WATTS 

NO 

- 

3 


XQ. 

= 

6.00000 E +01 

CM 

FQ 

CZ 

-l. oao ocE+on 


xo 

— 

6 . 0 0 0 0 u L + 0 1 

CM 


— 

0 . 


XQ 

H 

6.G0DG0E + ;U 

CM 

FQ 

= 

1.G3Q00E+GC 


NELEV 

- 

2 ' 


XELEV 

— 

0. 

CM 

ELEV 

= 

0. 

CM 



elevation point no, 2 

DISTANCE TO POINT XELcV = 

ELEVATION OF POINT...,*.. ELEV = 

INTEGRATION STEP SIZE OX - 


1 . 8G0CCL+ 02 
2* QQO C0E+ 3G 

2. 00QG0E+0Q 


ANOTHER CASE CC=NO, 1-YES1.... 

USE LAST POROSITY DISTN. ...... 

REQUIRE St LF 


NCASE 

L ASTPHI 
ISLFPRM 


j 

C 

0 



FINAL SOLUTION 


ro 

Co 


the maximum heat-tr^ns^er rate 

THE TOTAL LIQUID IN rflCK IS 


DISTANCE 

ST RE S S 

(CM) 

(CM LIQ.) 


0. 

1.1280 E+GU 

2 • 3) C u F + n o 

1. 1526E + 3C 

4.0330E+00 

1.17415+00 

6.0000F+G0 

1* 13845+G3 

3. QOCOE+OO 

1. 22365+00 

1.G0C3E+G1 

1*2 49 55+00 

1 .2UQ0F+ai 

1. 27635+00 

1.4QD0L+C1 

1*33405+33 

i *6000 c+Ci 

1. 332*45+30 

l.aoooE+oi 

1* 3617 5 ♦ G 0 

2 • 0 3 30E + G1 

1 . 39195+00 

2.200CE+Q1 

i . 42285+ 00 

2*4003 E +01 

1.4 5475+00 

2.6Q00E+01 

1.48735+30 

2.30G3F+U1 

1*52085+33 

3.03CGE+G1 

1 . 55515+00 

3. 2000E+01 

1. 59025+00 

3.43005+01 

1*62625+00 

3.6J0GE+G1 

1* €630 5+00 

3.S3G0E+01 

1* 7006E + U0 

4. 00C0E+G1 

1 • 73915+30 

4.23 COt +0 1 

1 . 77845+00 

4*43005+01 

1.81355+00 

4. 63005+01 

1 . 85955+00 

4. 100 35+ Oi 

1 . 9J 135 + 0 Q 

5.03005+01 

1.94395+00 

5.23 joE+01 

1 . 9 3 7 4 £ + 0 0 

5*4)035+01 

2 • 0 31 75+ w 0 

5*60005+01 

2. 37635+Cu 

5. 33005+01 

2. 12235+30 


S 1.154305+02 

WATTS 

5610 0E + 01 GRAMS 


STATIC HEAD 

PQROSI TY 

(CM LIQ.) 



1.12805+00 

8. 72 65E-G 1 

L. 1532E + G0 

a. 72 655- Q 1 

1 • 1 7 ?4 E + 0 0 

4. 72 65 5-01 

1.19475+00 

8. 72552-0 1 

1. 2169E+C0 

8. 7265 E -01 

1.2391E+00 

8 . 72o 5c- 3 1 

1.261 3 t+CO 

8. 7265E-C1 

1 • 2 8 36 E + Q G 

8. 7265E-G1 

1 • 30 58 1 + 0 0 

8. 7265c- Cl 

1.3230 5+0 0 

8. 72555-01 

1*15325+00 

8. 7265E-Q1 

1.3724F+C0 

8.72 65 5- 0 1 

1. 3947E+00 

8. 72665-01 

1.416 9 E+GG 

8. 72 655- C 1 

1.4391L+C0 

8. 72655-01 

1.46135+ 00 

8. 7 2655-01 

1.48 36 £ + 00 

H.'7265c-Ql 

1 .50525+00 

8 • 72655-3 1 

1.523C5+CQ 

3. 72665-01 

i * 5502E+G 0 

8. ?265c-01 

1.5 724 E +0 ii 

8. 7 2 6 5 5 - C 1 

1.5947E+00 

8 . 7 2 6 '> c - u 1 

1. 6169E+0Q 

8. 7265E-C1 

1.63915+00 

8* 7265c. - il 

1.66185+00 

•H . 7 2 6 5 E - 0 1 

i . 6 8 36 E + 0 0 

8. 7 2 6 5 E - 0 1 

1.70585+00 

8. 72655-01 

1. 72805+CO 

3. 7265 £ -01 

1.75025+00 

8. 72655-01 

1. 77245+00 

8. 7265E-01 


GRIT. STRESS 
(CM LIQ. ) 


2.2052E+OG 
2.2G52L+00 
2.2052F+3G 
2, 2G52E+0 0 
2.2G52E+0Q 
2.2052000 
2.2052F+G0 
2.2D52F+0C 
2.2C52E+0C 
2. 2 Q 52£ + 0 0' 
2. 2G52E + 0 0 
2. 2352 £4-00 
2.2G52E+D0 
2.2u52£+3G 

2.2G52r+QH 
2. 2G52E+ G C 
2. 2052E+GC 
2. 2 C 5?E + 0 0 
2.23520+0 0 
2. 2c. 52 l + 0 0 
2.2,: 52E+30 
2. 2u52F+ 3 G 
2.2352E+QC 
2.2352E+30 
2.2052c+ 0 0 
2 . ? C 52r + 3 D 
2. 2 0523" +0 0 
2. 20523+ JO 
2. 2 -525+ J 0 
2.2^522+00 



ro 

to 


6.0 10 J'c + 4l 
6.200 Jc+Gl 
6.4JoOE+Ql 
6.61 C 06 ♦ u 1 
6.8JC0E+Q1 
7 . 00005+01 
7 .20005+01 
7,40OCfc+0i 
7.60uGE+Cl 
7 . 3 0 0 0 E + 0 1 
ft, 3 1006 +G 1 
8.2JQ0E+O1 

0. 40005+01 

3.63035+01 
8. 80005+01 
9. GQC 36 + 01 
9 • Z Q 0 G E + C 1 
9.4100 £+01 
9.530GE+G1 
9*aocoe+oi 
1.00GQE+Q2 
1* 0 2 G 0 5 + 0 2 

1 , Q 4OOE + 0 2 
1* Q6C05 *-0 2 
l.GJGUt+02 
1.13005+02 
l.l2CuE+Q2 
1 . 140UF+O2 
16160 0E +02 
1*15 0 0 £ + 0 2 
1,23005+02 
1.220 15+02 
1.24005+02 
1, 26 DOE *0 2 
1.23005+02 
1.31 ODE *02 
1.32C0E+02 
1.34CGF+C2 
1.360 15 + 02 
1 • 3 3 G 1 E * u 2 
1*4100? *02 
i . 4 2 0 0 £. + C 2 


2.1696E+0U 
2 • 2 1 67E+ 3 3 
2*2 6 452* 3 U 
2. 3132E+00 

2 * 36295+00 
2 • 41 37E + GO 
2 • 465 6E+ 0 0 
2.5187E+30 
2*57295+00 
2. 62342*03 
2. 6ft52E*00 

2. 74335+00 
2.BS29E+09 
2. 8 6 4 0 E * 0 0 
2 • 92665+0 0 
2.99095*30 

3 • u 5 7 0 E * 1 3 
3. 12435+00 
1* 1946E+G1 
3*2664 £ *00 
3.3. 40 45 + 03 
3*41b6E+0O 

3 • 4 95 2 E* 0 0 
3* 57635+ 03 
3*66025*00 
3.7463E+30 
3 . 83665+00 
3.92955*00 
4 • G 25 95 + 0 0 
4, 12605*00 
*♦. 23015*03 
4.33745*00 
4.44615*03 

4 * 5 5 61 E+ o 3 
4*6672E+C0 
4,77935*00 
*4.59235 + 00 
5 . 10535+00 
5*11975+03 
5.23195+00 
5. 34795+03 
5.46155*00 


1.79475+00 
1. 31695+00 
1. 33915+Cu 
1*8613 5*00 
1* 88365 + 0 0 
1.90585+00 
1. 928G6+CG 
1 . 9502E + 0 0 
1.97245+CG 
1.99475*00 
2.01695*00 

2.0 39.1 E + GG 
2.06135*00 
2.08365+00 
2. 10 5b £*00 
2.1233E+00 
2.1502E+C0 
2.17245+OQ 
2.19475+CC 
2. 2169E+0D 
2.23915*30 
2. 20135+00 
2.283.6 5 + 00 
7*305 8 E+GG 
?.3230f+00 
2. 55U2E+0C 
2. 37245+00 
2* 3 947 E +00 
2.41695*00 
2.4391E+0O 
2.4613E+Q0 
2.48365+00 
2. 50585+00 
2. 52305+00 
2. 55O2E+G0 
2. 57245+00 
2. 59475+CO 
2. 61695+OC 
2. 63 91 5 + 00 
2.66135+CG 
2*68365+00 
2.70585+00 


8. 7265E-01 
8. 72075-01 
8.6967c.— 01 
8. 67245-C1 
8. 6477E-G1 
8 *62275-01 
8. 5972E-01 
8. 57135-01 
8. 5453E-S1 
8. 5183E-01 
8, 49H5-.0.1 
9, 4635E-01 
6 .4353E-C1 
8.40675-01 
9. 37755-01 
9.34775-01 
8. 31745-C1 
8. 28 645-C1 
8. 2543E-01 
8. 222&5-G1 
8. 1396c. - G1 
8. 15 595- Li 
e. i?i55- r i 
o. 08625-01 
ft. 05015-01 
8. 01315-11 
7. 17525-01 
7. 13625-01 
7. 89625-01 
7. 3552E-C1 

7. 9129- -01 
7, 76995-C1 
7. 72665-01 
7. 6 9 3 45 - G 1 
7. 6412 5 -Cl 
7. 69715-01 
7.554' 35- u 1 
7. 51165 -Cl 
7. 46935-01 
7 . 4 2 7 4 K - G 1 
7.38635-11 
7 • 3 4 5 2 E - C 1 


2.205 25+0 0 
2.2167E + G3 
2.26456 + D C 
2.3132^+QQ 
2. 3629E+0 C 
2. 4137E + 0 0 
2.4656E+Q0 
2.51875+00 
2.57295+10 
2.628 4£ +00 
2.6852E+Q0 
2,7 433E + 00 
2.8 1295+ G 0 
2. 364CE+C 0 
2.92665+00 
2.99095*00 
3.35705+00 
3.12485+00 
3. 19465+00 

3. 2 6645+ 3 C 
1.34 n 4F+0r 
3. 4 166E+ 1 C 
3.495 2 E+30 
3.57636+0° 
3.66025+00 
3.74685 + 0 0 
3.8366E +0T 
3.929 5F + 0 0 
4.02595+30 
4.1260^+00 
4.23315*30 
4.33746+00 
4.4461 E+Q H 
4* 56618+0 0 
4.66725+00 
4.7793E+Q0 
4.69235+00 

5. 0 05 8E + 0 0 
5. 11975+00 

5.2 339F + CC 
5.3479E+0G 
5.4615E+CG 



1.44G0E+G2 
1,4500 E. + 02 
l*4*00t+02 
1.50G0E+02 
1.52GGE+02 
1.54007+02 
l*56u f 3E + G2 
1.53C0E+02 
i.60 0 3 E + 0 c 
1.620GE+Q2 
1 » 64 0 0 E + 0 2 
1.66CGE+G2 
1 * 63 u 0E + G 2 
1.7300 E +32 
1 • 720 3 £ + 0 2 
1*740 QE +02 
1 .760QE+02 
1 *73 0 0E + 0 2 
1*3300 E+ 02 


5.5744E+G3 
6 • 6862E+ 3 3 
5* 7966E+G0 
5 •9B52E+GQ 
6 • Oil 5 E+ 0 3 
6. 115GE+00 
6*21 54E + 0 0 
5.3122E+00 
6 • 40 4 3E+ J 3 
6 • 492 3E + 0 3 
6*5 75 6E+ u 0 
6*652 3 E+03 
5*7233 E+ 00 
6 • 7 634 1 + 0 0 
b • 845 9£ + 0 3 
6* 8961E+u9 
b • 933 7£+ 0 3 
6*9732 £+00 
6*9996 E+Oj 


2.72'JOE + OQ 
2. 7502E+GC 
2. 7724E+E0 
2.7947E+CQ 
2*3 169E + P 0 
2. 3391E+GG 
2* 3613E + G0 
2* 3636E + 00 
2.9058E+03 
2 • 9 ? 3 u E + £ G 
2.9502E+GG 
2.9724E+0Q 
2* 9947E + 0G 
3.3169E+00 
3.0 391E + CG 
3. J616E+GQ 
3.Q336E+Q0 
3.1G58E + 00 
3.1233E+CG 


7* 3Q52E-01 
7 • 2659E-G1 
7. 2275E-C 1 
7. 19Q2E-C1 
7. 154GE-01 
7. 1191E-Q1 
7*0 3561: - 01 
7. 0 536E-01 
7. j 2 32E-C 1 
6*99465-01 
6. 9679E-C1 
6. 9432E-01 
b. 92Q6E-U1 
6* 90Q2E-01 
b* 3 321E- 0 1 
6 • 9 6 b 4 E - 0 1 
b. 3 5 3 2 £ - C 1 
6*3 424E - D 1 
6. 3 34 3E- u 1 


5*5 744E + 09 
5* 6862E + 0 C 
5* 7966E + U 0 
5.9G52E+30 
6.011 5E+0C 
6. 1150E+00 
6 • 2154E+ Q n 
6.3122E+QC 
6*404 3E +00 
6 • 4 9? 3E + 9 0 
6.5756E+3C 
6*652 3E +30 
6*723 3 E+GQ 
6.7 384E + 0 C 
6* 3459E + 0 0 
6.8961E+Q Q 
6.9337E+00 
6.9732F+0C 
6. 9996E+0C 



GRADE SAMPLE PROBLEM CASE NO. 4 

NON* SELF- PRIMING WICK ZERO GRAVITY 


INPUT VARIABLES AND FLUID PROPERTIES: 


LIQUID NUMBER 

TEMPERATURE 

LIQUID DENSITY 

SURFACE TENSION 

LIQUID VISCOSITY 

LATENT HEAT.. 

GRAVITATIONAL ACCELERATION. .. . 

WICK AREA 

WICK HEIGHT 

WICK p I BER DIAMETER 

INITIAL POROSITY 

FINAL POROSITY 

FINAL STRESS 

VAPOR-SPACE DIAMETER 

NOMINAL HEAT-TRANSFER RATE.... 

NO. HEAT-INPUT SECTIONS 

SECTION NUMBER 1 

SECTION LENGTH 

HEAT-INPUT FRACTION...... 

SECTION NUM3ER 2 

SECTION LfcNGTH... 

HEAT-INPUT FRACTION 

SECTION NUMBER I 

SECTION LENGTH. . , . 

HEAT-INPUT FRACTION 

NO. ELEVATION POINTS 

ELEVATION POINT NO. 1 

DISTANCE TO POINT 

ELEVATION POINT 


LIG 

= 

2 

TKELVN 

= 

2 • 943 C C E ♦ G 2 

RHOL 

= 

6.C9A35F+C2 

3T 

= 

2.094G9E-G? 

VI SL 


1.37S4CE-Q4 

HED- 

— 

1.13576E* 0 6 

GES 

= 

0 . 

AW 

- 

4.C1QUGL-31 

HW 

- 

1 • 123 t iiE + 0 0 

DI A r 

- 

1.C16Q0E-Q2 

PHIQ 

= 

0 . 

PH IF 


0 . 

PF 

• = 

7.Q303GE+G3 

QIAVS 

= 

4. 7-3 3 u 0 E- 0 1 

ROOT 

~ 

1 . GG30GE+Q2 

NO 

— 

3 

XQ 

Z 

6.Qa0C0E+Gl 

FG 

~ 

-1. QOOCCFfDO 

XG 

Z 

6 • 0 3 0 0 1 L + D 1 

FO 


0. 

XQ 

z 

6.0000 CE+01 

FQ 


1.0000 GL+QO 

NELEV 

= 

2 

XELEV 

z 

f . 

ELEV 

- 

0. 


DECREES KELVIN 
KG/CU. M 
NEWTONS/M 
NEWTON-SEC/SR. M 
JOULES/KG 

standard gravitt 

SQ. CM 

CM 

CM 


CM LIQ.. 
CM 

WATTS 


CM 


CM 


CM 


CM 

CM 



ELEVATION POINT NO. 2 

OTSTANCt TO POINT XELEV = 

ELEVATION OF POINT ELEV - 

INTEGRATION STEP SIZE OX = 

ANOTHER OASE It- NO, 1=YESI.... NCASE = 

USfc LAST POROSITY DISTN....... LASTPHI = 

REQUIRE SELF PRIME • ISlFPRM = 


!,SOaOQE4-u2 Oil 
a.iioaocE + oo cm 

2.QOOOCE+OC CM 


1 

1 

n 

c 



FINAL SOLUTION 


THE MAXIMUM HEAT-T RANSFER RATE IS 2.0QijOOF402 
THE TOTAL LIQUID IN NIC< IS J.56100E*0l 0RAH3 


LO 

CO 


DISTANCE 

COM) 


0. 

2,ooooe*oo 

4.0i3Uut>ijO 

6.000CE400 

3. (JGUu£4Q0 
I.0fl00E»01 
1.2000E+0I 
i.40G0£+Cl 
1 • 6 ] 00E ♦ C 1 
1.3JD0E+C1 

2 . 0 0 0 C E 60 1 
2* 23uuE+Gi 
2 • 40 U OF. *01 
2 . 6 QQPE 401 
2.fiQG0£ + 01 
3.00G0E4G1 
3.2G0CE+G1 
3 • 4000E+G1 
3 • 6 3 u uE + G 1 
3.8J0CE4C1 

4. JQOOEtGl 
4.200 uE4Ci 
4.43i;3F40i 
4.6U0CE+Q1 
4. 800 QE 401 

5 . 0 3 0 0 E 4 0 1 
S* 23 C 0 E*C 1 
5.43uoE+Cl 

5.6JJ0E4C1 

5.3000E+C1 


STRESS 
(CM LIQ.I 


2 • 9340E- 3 1 
2 • 9912E-0 i 
3* C129L-G1 
3.D490E-01 
3.Q995E-01 
3.1645E-01 
3.2439L-01 
3.3377E-01 
3* 4459 E- 01 
3.5646 E- 01 
3. 7053E-J1 
3. 3573E-Q1 
4. 0233 E-Gl 
4*203 3 E- 01 
4 • 393 7E- 0 t 
4*608 3E-0 1 
4. 8317E-01 
5 • G699E-0 1 
5 • 3225 F- 01 
5*5395 E-Q i 
5* B71-JE-01 
6* 1669F-01 
5 • 4773 E-Ul 
6 . 80212-0 1 
/. 1413E-G1 
7. 4950E-C 1 
7 • 66 31c.-* i 1 
3 « 2 4 5 6 F - 1 : 1 
3 • 64256-01 
9*053 9 E-01 


STATIC HEAD 
(CM LIU* > 


0 . 
0 * 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 * 
0 * 
0 . 
0 . 
3 . 
0 . 
0 . 
3* 
0 . 
■J t 
0 * 
0 . 
1 . 
0 . 
0 . 
a. 

o . 


WATTS 


POROSITY CRIT . STRESS 

<CH LIQ. ) 


3. 7265E-01 

2.2052E40G 

3. 726 5 E- Cl 

2.2052E+00 

6* 7265E- C 1 

2.2 052E 4 Q 0 

8. 7265 E -Cl 

* 2. 2052E+0G 

3. 7265E-G1 

2. 2 0 525.4 0 G 

8. 726SE-01 

2.2052E4QC 

8. 72655-11 

2.2f!52E400 

3* 7265. E - Cl. 

2 . 2 C 52£4 0 G 

3* 7265E-31 

2.2052E+0G 

0* 7265E-G1 

2*20522400 

3 • 72 65 E- 01 

2. 2 0 52£ 4 Q G 

S. 7265E-Q1 

2*205 2h 400 

8* 7265E-U1 

?#2u52t+0C 

8. 7265E-01 

2. 2u52£4Q G 

8. 7265E-01 

2. 20526400 

8*72 652- L 1 

2.2J52L+GC 

3. 72652-ul 

2.20526400 

8. 72 652-01 

2*20 52£ 4 G 0 

8. 7265E-D1 

2.2052E40G 

8* 7265E-Q1 

2*2-52 £400 

8* 72 65 2- L 1 

2. 20 52 £4 QC 

6* 7265L-Q1 

2*2 i :52f.:40C 

8. 7265F-C1 

2.2u52E40n 

8. 72652-01 

2.2052640 G 

8. 72652-01 

2.2G52E+CC 

8. 72y5E-ri 

2 • 2 C 5 ? l 4 0 0 

3. 7265E-C1 

2.2-52E40C 

a * 7265n — ul 

2.2052E40 0 

6* 72o5E-01 

2.2 052^40 0 

8* 726 52-01 

2.265?£4j0 



6 • 0 3 0 QC + C 1 

9 .479SE-G1 

6.2COOE+C1 

9*9124 E- 01 

6.48Q0E+G1 

1.03552+03 

6.630GE+C1 

1.0315E+00 

6.80002+01 

1.1292E+00 

7.0 JCQE+01 

1 • 1787E + 0 0 

7*2300 E+ Cl 

1 • 2 33 It + 0 0 

7.4QC0E+01 

1*233 52 + 0 0 

7 .bOGuE+ul 

1 • 3390E + C 9 

7.83C0E+ 01 

1*31662+09 

d.flQ.OOE + Ol 

1 .45652+00 

8 . 20002+01 

1 * 5187E+ G 9 

8.40CGE+G1 

1. 5335E+00 

S.69Q0E+C1 

1 • 6 5 u 3E + 0 0 

8.80002+01 

1* 72392+09 

3*0JiiCE + G 1 

1. 7 9 5 3 2+00 

9 • 23002+01 

1* 86972+33 

9 • 4 JO 3E + Qi 

1.94332+03 

9.63302+01 

2.0 31 22+310 

9.3J&32+C1 

2.11712+33 

1.09G0F *02 

2.236 72+ 0 J 

1.0 2002+02 

2.30332+33 

1* 04002+02 

2.3933E+G0 

1 . 06 Out + C 2 

2 • 5 3 3 1E+ U J 

1.O803E+C? 

2. 69632+30 

1.13002+52 

2 . 71352+09 

1 • 1 2 C J F + C 2 

2 . 8 35/52+ 0 J 

1 . 1 4 G 3 E + C 2 

2.9583E+LJ 

t.i60QE+G2 

3.03662+39 

1.13G3E+G2 

3.22152+33 

1. 200GE+G2 

3 . 36542+03 

1.2200F.+C2 

3*51382+ 0J 

1.24QQE+02 

3 . 66072+03 

1 • 2 6 0 0 2 f G’2 

3. 81272+00 

1.28012+02 

3*96682+03 

1.33GG2+02 

4.12252+03 

1. 32DQC+Q? 

4 . 2 7962+ 9 3 

1 • 3 4 G 0 F + 02 

4.4379E+3J 

1.3600E+02 

4 . 5 1 6 3 h + G J 

1.3300 £+02 

4 • 75602+09 

1.43032+02 

4.91502+03 

1.420UE+C2 

5. 0 75 42+03 


o o j ij o o o cc □ "3 c o c □ q o ti ca o ij o q o □ a c o o c a o o a c o j o 


0. 

8 • 72 65c.* 11 ' 0 1 

2.2G52E+GG 

9. 

3. 72972-01 

2. 21672+30 

0. 

8.69672-01 

2.2645E+0C 

0. 

8.6 724E-G 1 

2.31322+GC 

a. 

8. 64772-01 

2 • 36292 + 0 0 

a . . 

8. 6227E-R1 

2. 41372+30 

a. 

8.5 9 722-C i 

2.46562+90 

0. 

8.57 1 3E-G i 

2 . 51372 + G C 

a. 

3. 54502-01 

2*57292+00 

0. 

G • 5 1 83E-G 1 

2.6284E+0C 

0. 

• 6. 49112-01 

2*63522+90 

0. 

8. 46352-01 

2.7433E + G 0 

0. 

8* 4353E-01 

2.8329E+0C 

a* 

8.40672-01 

2.864 Dt+GC 

a . 

ft. 37752-G1 

2. 92662 + 0 0 

a. 

ft. 34772-01 

2.99092+00 

0. 

8. 31742-01 

3.057 0E +00 

3 * 

ft. 28642-01 

3.1 2482 + G C 

0 * 

ft. 25482-01 

3.194 6E +00 

Sf * 

8.22262-01 

3. 26642+ 0G 

3. 

ft. 18962-01 

3. 34042 + 0 G 

G. 

tt . 1554.2-0 1 

3.4166L+G0 

0. 

8.12152-01. 

3.4952E+0C 

0. 

ft. nft622-Ci 

3.5763E + J C 

u * 

8.G5012-C1 

3 .66022 + 0 C 

0. ■ 

3 * 31312-31 

3, 74632+0 D 

u. 

7. 97622-01 

3.83662 + u 0 

n. 

7. 3362F-C1 

3.9295E+C f« 

0. 

7. 8962E-S1 

4.97592 + 0 0 

0 • 

7. 85522-C1 

4.12602+00 

0 ♦ 

7. 81292-*! 

4.2301E+0G 

3 . 

7. 76 i 8F-G1 

4.3 3745+'"' 0 

a. 

7 • 72662-01 

4. 44612 + U 0. 

o • 

7*68342-01 

4. 55612+00 

3 . 

7. 64J2 z- r i 

4.56722 + ^ 

0 . 

7.59712-01 

4. 779 32 + DC: 

r* 

v • 

7, 55432-01 

4.8923E+C0 

U ft 

7.51162-01 

5.0 05 82 + 0 C 

-J. 

7. 46932-01 

5.11972+00 

' J ft 

7. 4 2 742-01 

5*253 9 r+C r 

9. 

7 * 3<5oD— . — cl 

5.3474P + j 0 

ft 

7 • 54572—01 

5.4 6152 ♦ 0 0 



1*4400E*Q2 

5.2135E+33 

0* 

7.3052E-C1 

5* 5744E* 0 0 

1.4600E+02 

5.3857E+00 

0* 

.7* 2659E-01 

5* 6362E + Q 0 

1.430QE+02 

5* 5385E + 00 

0. 

7*22 752-01 

5.7966E+00 

1. 500GE+G2 

5 * 68 3 QE* 00 

0* 

7.1902E-O1 

5.9052E+00 

1 • 5 2 0 0 E * b 2 

5.3337£*30 

3* 

7 * 1 540E-0 1 

6* ti 11 5E + 0 0 

1.540QE+G2 

5.9746E+GJ 

a* 

7*1 1 91E-01 

6*1150E+QC 

1.56GOE+02 

6.1101E+30 

0. 

7*0856 E -01 

6 • 21 54E + 0 0 

1.5300E*02 

6.2392E*G3 

0 • 

7*0 5 36E-0 1 

6* 3122E + 0 0 

1*600 OE +02 

6.3611E*0a 

0* 

7.Q232E-01 

6.4048E+Q0 

1 *6 20QE + G 2 

6* 475 0E*O 3 

0* 

6, 9946E-C1 

6* 4928E. + 0 0 

1.6400E+G2 

6.5800E+C3 

3* 

6* 96 79E-0 1 

6* 5 756E + G0 

1 • 6600E + 0 2 

6* 6752 E +00 

0* . 

6. 9432E-01 

6*652 8E +00 

1*680 QE* G 2 

6*75986+00 

0* 

6 * 92C6E-0 1 

6.7238E+00 

1*700 QE +02 

6*8 3 31F* 0 0 

0 . 

6* 90C2E-01 

6.7834E+0C 

1.720QE+02 

6.8943E+03 

0* 

6. 3 82 IE- 01 

6.8459E+C0 

1.74UQE+02 

6 • 9428E+33 

0* 

6. 3664E-L1 

b* 3961E + 3 C 

1.7600E+02 

6 • 973 Q E + G G 

3. 

6. 3532E-01 

6* 938 7E + C P 

1.730QE+G2 

6 *'9993E+ 30 

0* 

6. 3424E-C1 

6* 9732E + 0 0 

l.aOOQE+02 

7 * 0 3 6 5 E* 0 0 

0. 

6. 8343E-C1 

6 . 9996E+ 0 G 



GRAOE SAMPLE PROBLEM 
HOMOGENEOUS WIC< 2-CM ELEVATION 


INPUT VARIABLES AND FLUIO PROPERTIES* 


LIQUID NUMBER**** 

TEMPERATURE * * 

LIQUIO DENSITY 

SURFACE TENSION 

LIQUIO VISCOSITY 

LATENT HEAT 

GRAVITATIONAL ACCELERATION.... 

NICK AREA 

WICK HEIGHT 

WICK FIBER DIAMETER 

INITIAL 3 GROSI TV 

FINAL POROSITY 

to 

cn 

FINAL STRESS.*. 

VAPOR-SPACE DIAMETER 

NOMINAL HEAT-TRANSFER RATE.... 

NO. HEAT-INPUT SECTIONS 

SECTION NUMBER i 

SECTION LENGTH 

HEAT-INPUT FRACTION 

SECTION NUMBER 2 

SECTION LENGTH 

HEAT-INPUT FRACTION. . . . . . 
SECTION NUMBER 3 

SECTION LENGTH 

HEAT- INPUT FRACTION 

NO. ELEVATION POINTS..... 

elevation point NO. 1 

DISTANCE TO POINT 

ELEVATION OF POINT 


CASE NO. 5 


LIQ 

- 

2 


TKELVN 

= 

2.9400CE+02 

DEGREES KELVIN 

RHOL 

- 

6.09475C+02 

KC-/CU. M 

ST 

= 

2. Q94C9E-Q2 

N EWTON S/M 

VISL* 

— 

1.37340E-Q4 

NEWT0N-SEC/5Q. M 

HFG 

- 

1. 1S576E+06 

JOULES /KG 

GEE 

= 

i.acooctMio 

standard graviti 

AW 

= 

4.G1Q00E-01 

3Q. CM 

HW 

= 

1 . 123 G CE* QC 

CM 

OIAF 

- 

i.Gl&DGb-O? 

CM 

PHIG 

- 

6.8350GE-3I 


PH I F 


0. 


PF 

- 

c. 

CM LIQ. 

DIAVS 

- 

4.7QO0OE-C1 

CM 

QOOT 

ZL 

l.QbflOCE+3? 

HATTS 

NO 

= 

3 


XQ 

- 

6*G00'30t*Ji 

CM 

FQ 

= 

-l. oooccr + ac 


XQ 

- 

e . c a o o cr > c i 

CM 

FQ 

= 

0 . 


XQ 

— 

6. OG 0 0 Gb 4- 01 

CM 

FQ 

= 

1 • 0 3 0 C 0 E + 3 0 


MCLEV 

r 

2 


XEIEV 

- 

0. 

CM 

EL EV 

- 

0 . 

CM 



ELEVATION POINT NO. 2 

DISTANCE TO POINT. XELEV = 

ELEVATION OF POINT ELEV = 

INTEGRATION STEP SIZE.... OX = 

ANOTHER CASE <o=NO, t=YES».... NCAS? = 

USE LAST POROSITY DISTN. ...... LASTPHI = 

REQUIRE SELF PRIME ISLFPRM = 


t*> 

-4 


1.8GQCC.E+G2 CM 
?. OOOOOE+OD CM 

E.QGDOCE+uD CM 


C 

0 

0 



FINAL SOLUTION 


THE MAXIMUM HE AT -TRANSFER RATE IS 2.93828E+01 WATTS 
THE TOTAL LIQUID IN WICK IS 3.30665E+G1 GRAMS 


<X> 


OISTANCE 

(CM) 


а. 

2.0D00E4-CG 
<»• dd'QGE + OO 

б. 03G0E+0U 
S.OOGQE+OO 
l.OQCOt+ul 
1.23Q0E+01 

i.Aaoae+ci 

1 ,63nn£fOi 
i .33002+01 
2 • 0 3 09E+0 1 
2.2QG0E+01 
2.400 OE+Oi 
2.63Q0E+C1 
2. 30032+01 
3.Q300E+Q1 
3.20002+01 
3.4300E+01 
3.60QGE+G1 
3.800GE+01 
4.0GG0E+01 
4.23O0E+OI 
4.4QGGE+0T 
4.6300 E + 01 
4.3 J002+Q1 
5.0000 E +01 
5.2GC0E+G1 
5.4GC0E+C1 
5.6000E+01 
5. 833QE + 0 1 


STRESS 

(CM LIO.) 


1.1280E+Q0 
1. 1 5l 3 E+ u 0 
1.17682+00 
1. 20442+00 
1.2342E+0S] 
1.26612+03 
1.30G22+00 
1 .3365t+GQ 
1. 37492+00 
i. 41542+00 

1 * 45S2E+ 00 
1.50312+03 
1 . 550 12+ GO 
i • 5993E+ G G 
1*65372+90 
1.7342E+G9 
1 .75992*00 
1 . 61776+03 
1.3773E+00 
1.939 9 E+30 
2*u342E+30 
2.0707E+03 
2.13942+00 
2 . 2102 2 +00 
2*2331 E+ 00 
2. 35322+00 

2 • 43552+ 0 3 
2.5149 E+uO 
2 • 5965E+ 0 0 
2. 6803E+00 


STATIC HEAD 
(CM LI Q • ) 


1.123GE+0G 
1.1532E+0O 
1. 1724E + 00 
1. 19472*00 
1 • 2169E+0G 
1 » 2 391E> GO 
1. 2613E+0G 
1.2836E+CG 
i. 3Q562+0U 
1.329'0E+00 
1.3502E+C9 
1.3724E+Q0 
1* 3947E + Go 
1.4169E+GQ 
1.4 39 1 E + G O 
1. 4613E+00 
1 . 4 ft 36 1 + G G 
1.50 5i3fc*0G 
1.5233E+C0 
1.5532EMJG 
1 .5724E+CQ 
1.59472+OC 
1 • 6169E + G C 
1.6391E+0O 
1 . 661 j»£ + 0 0 
1 . 6 8 36 2 + G 0 
1.7C58E+U0 
1. 72 3CE+GC 
1. 7502E+OO 
1. 7724E+G 0 


P9R0SI TT 


6. 83532-01 
6. 83506-01 

6. 33502-01 

6.83502-01 

6. 83502- 01 

6. 33502-01 

6.83502- d 
b. 38502-01 
6* 83 5 3E- u 1 
o. 63502-01 

6. 33502-01 
6* 83502-01 

6. 33502-01 
o*83 502- w 1 

6. 33502- 01 
6. 8353E-U1 
6. 8350E-O1 
6* 3 3502-01 

6.83502- 01 
6 • 83502-01 

6. 33502 - Cl 
6 . 8 3 5 0 n — v 1 
6. 8 3 b 0 i — 01 

6.33502- 01 
6. 335o2-Ul 
6 • 33502-01 
o * 8 3 6 0 2-0 1 
6 .635 0 2-01 
6. 3 3 5 J 2 - 0 1 
6 * 33502-01 


CRIT • STRESS 
(CM LI Q.J 


6.9972E+Q0 
6 . 9 9 72 £ + 1 G 
6.9972E+00 
6. 99722+00 
6. 9972.2 + 0 0 
6 • 9972E + 0 0 
6. 9972E+0 0 
6.9972E+30 
6. 9972E+0C 
6. 9972E+H0 
6.9972E+QP 
6. 9 972E + G C 
b. 99722+ 0 C 
6. 9972F + G 0 
6 . 99722+80 
6. 9972E+0 G 
6.9972E+3G 
6.9972E+00 
6. 99^22+0 0 
6.9972E+0G 
6.99722+00 
6 . 99 7 2E + 0 C 
6.9972E+00 
6.99722+00 
6.9972 E+OC 
b. 9972- + 0 C 
6. 99722 +0 r 
6. 99722 + 3 G 
6. 99722 + 3 C 
6.997PE+0C 



u> 

KO 


6.0003^+01 
6.2JGJE+G1 
6 . 4 3 G 0E + 0 1 

6 • 6 3 0 ii £ * 0 1 
6.aaoaE+oi 
7.00006+01 
7.2300t+0i 
7.49C3E+C1 
7.6000E+01 

7 * 8QQDE + G 1 
8 • Q 0 GuE + O 1 
3. 2 D0GE+01 
8.430GE+31 
8.60006+01 
B.800GE+01 
9* Q0&GE+Q1 
9.2OGOE+01 
9.40aOfc+Gl 
9.60G3E+G1 
9.83uGE+0i 
1.GOGOE+D2 
1*021 OE+G 2 
1 . 340GE + 02 
1 • 06 DOE + G 2 
i.o joue+02 

1.130G6+C2 
1*120 OE *02 
1*14GQE+U2 
1 .160 Ot ♦ U 2 
1 • 1SGGE + G 2 
1.23QGE+G2 
1.2203E + G2 
1. 24006 + 02 

i.26C0E*“Q2 
1.2800 E+G2 
1 .3300E+G2 
1*320 GE+32 
1.3 4G0F+C2 
i# 36G0E+G2 
1.39G0E+02 
1.43002+02 
1.42CQE+G2 


2. 7662C+Q0 
2#8523l+ 0G 
2. 93956+30 
3. 0261E+JQ 
3* 1127 E *-33 
3. 1994E+GD 
3.2360E+uG 
3. 3726E+Q J 
3 • 45922+ 0 3 
3 . 5459E+ uu 
3 • 6 325E + G Q 
3 • 7191E+03 
3. 60546+03 
3.892 4 E+ CD 
3 .979JE+0J 
4.86576+03 
4 . 1523E+0 J 
4.2389E+G0 
4. 32566+ 0 0 
4 . 4122E + 0 D 
4.4 98 8E + 3 3 

4* 53556+03 

+ *67216 +uG 
4. 733 70 CO 

4 • 84536+GO 
4.932DE+C3 

5 . C136E+C3 
5*iJ52F+GD 
5 . 1919E+G J 
5.2785E+G0 
5. 3551E+C J 
5.451 JE+OQ 
5. 5343E+G3 
5 • b 16 4E+ 0 3 
5.6958E+G0 
5 • 77 31E+ u 3 
5 . 8482E+03 
5 • 9212E + u 0 
5.992QE+UG 
5 . 3 6 0 6 E + 3 0 
6.12710 0 0 
6 .19146+30 


1 • 7947E+ 0 0 
1.81696+00 
1* 8391E + GG 
1*86136+00 
1 • 3 8 36 E ♦ G G 
1. 90536+00 
1.928GE+00 

1 * 95 3 2E + 0C~ 
1.9724E+0G 
1.9947E+00 
2* G 1696 * C 0 
2.0 391E ♦GO 
2.06136+0 0 
2.G336E+00 
2.10 58E + u Q 
2.1230 £+03 
2.15G2E+CG 
2.1724E+0C 

2 • 1 9 47 E + Q0 
2 • 2169 1 + CO 
2*2391 f +00 
2* 2613 E +00 
2.2836E+GC 
2.3058E+00 
2. 3230E+0C 
2.35326+00 
2.37246+00 
2* 3947E + QG 
2.4169E+0G 
2.4391E+00 
2 • 46136 + 0 0 
2*48366+00 
2.5G58E+0Q 
2.523GE+Cu 
2. 55D26+CG 
2*5724 E+Lu 
2. 5947E+C0 
2*61696 *00 
2.63916+00 
2.6613E+0C 
2 * 6 83b£ + 0 0 
2. 70586+00 


6. 335GE-01 

6. 33506-01 

6.33506-01 

6. 83506- 01 

6. 33506- 01 
6. B35DE-01 
6* 835QZ-01 
6* 335 JE-ul 

6. 33506- C1 
6. 8550E-C1 
6.3353E-D1 
6. 3353E-01 
6* 335 j£- ul 
6* 335CE-01 
6* d?50£-Oi 
6. 335 JE-Qi 
6* 3350E-01 
6.8350 E -Cl 
6. 335GZ-C1 
6. 3350E-C1 
6. 3350E-C1 
6. 3 3 5 0 £ - 0 1 
o* 33535-01 
6* 8 35 0 E-fc 1 

6 . 83506- 81 
6 * 33506-01 
b. 83536-PI 
6.835Qc-Cl 
6. 3 3 5 3 E - 0 1 
6. 3 3 5 u £ - J 1 
6 * 3350E-01 
6. 8350E-C1 
6* 8 3 5 0 E - 0 1 

6.33506- 31 
6. 9350E-01 
6* b3 6QE-Jl 
6. 335 JE -01 
6* 8 3 5 u£ - u i 
6.8351c- J1 
6. 3353E- 41 
6. 835GE-C1 

6.33506- 01 


6* 9972E + 0 0 
6*99726+00 
6 • 99726+00 
6.99726+00 
b • 9972E + 0 C 
6*99726+00 
6.9972E+GC 
6.9972E+DP 
6 • 9 9 7 2 1 + 0 P 
6.9972E+H0 
6.9972E+Q0 
6.9 9726+00 
6.9972E+00 
6*9972 E+00 
6.99726+QG 
6.99726+00 
6*9 9726+0 0 
6. 9972E+CC 
6.99726+00 
6. 9972 r+C 0 
6. 9972 1. + 3 i 
6*9972 E+00 
6*99726+30 
6.997^6 + 3 C 
6*99726+00 
6 . 9972E+0C 
6.9 972E+0C 
6 . 9 9 7 2 E + 9 C 
6. 9972E+0 0 
6.91726+0° 
6. 9 3 7 2 1 + 9 0 
6.99726+0 [■ 
6. 9 J72E + C 0 
6*93726+00 
6*99726+00 
6*99726+30 
6.99726 + 0 C 
6. . 9 9 7 2 1 + H C 
6. 99726 + 0 0 
6.9972E+QC 
6.99726+GC 
6.9972E+O0 



1.440GE+O2 
1.46GJt+G2 
1.48Q0E+G2 
1.50002+02 
1.5200L+C 2 
i • 54 0 QE + 0 2 
1.56QOE+C2 
1.5300E+02 
1.69002+02 
1 • 6 2 G 0 2 + 0 2 
1.64C0E+02 
1 .66G9E+C2 
1.63 COE + 0 2 
1*7 j G 3 1 + G 2 
1 . 7 20 Q2 + G 2 
1.74002+02 
1. 76 0 Of* 02 
1.74 DOE *0 2 
1. 30002*02 


6* 25j6E*Q Ej 
6.31 36 E+ 09 
6 « 37142+ 0 0 
6. 4’7iE+Q0 
o*43 0 6 E-*QQ 
6.53202+30 
6, 53122+00 
6.62 33 £+ Q Q 
6.67322+00 
o. 71532+00 
6 • 75652+ 9u 
6.79492+ 30 
o.vHllt+JO 
6* 36522+0 3 
6 * 89 72 2+ 0 J 
6.92692+03 
o. 95462+00 
6.93002+00 
7.03332*00 


2 • 7239 E ♦ 0 0 
2.7502E+00 
2.7724E+00 
2. 79472+09 
2 • 81692 + 0 G 
2.8391E+0C 
2.3613E+0G 
2. 8B362+00 
2.9G582+CG 
2.9230E+00 
2.95022+00 
2. 9724E+D0 
2. 99472 f O 0 
3. 0169E+0G 
3 . 0 391 £ + 00 
3.36132+00 
3. 38 36 £+00 
3 • 1 G 56 E + G 0 
3. I239E+C0 


6. 33532-01 
•6. 335 02- oi 

6.33502- Cl 

6. 3 3502- G 1 
6 * 8 3 5 9 £— b 1 
6.3350E-G1 
6. 335Q2-C1 
6. 8350--01 
6. 8359 £-01 
6. 33592-01 

6. 33502- 91 

6.83502— 01 

6. 83502- 01 
6. 3359*1-01 
6* 8 3502-01 
6. 8 3 5 0 2 - Q 1 

6. 33502- 01 

6.33502-01 
b. 81502-01 


6.9972E+0 0 
6. 9 9 722 + Q C 
6. 9972E + 0 0 
6. 9972E+ 0 0 
6.997?£+00 
6. 99722+0 C 
6.9972E+39 
6.9972E+0G 
6.99722+0 fi 
6.9372E+00 
6.9972E+00 
6.9972E+3C 
6. 99726+09 
6. 9972 E + 9 0 
6. 99722 +00 
6.9972E+90 
6. 9972E + 0 G 
6.99722+90 
6 • 9972E + 30 



APPENDIX B 


PROGRAM LISTING 


41 



RUN* COMPILER (VER.2.3M) 


03 / 23 / 74 . 15 * 31 * 43 * 


G u8 004 
0QQ0C4 
0 G00G4 
DC30G4 
G GO 094 
003004 

Q 00004 
0C0GG6 
0 09 OU7 
G 0001Q 
0 00012 
000014 
000016 
8 00 020 
0 00 0 22 


PROGRAM M4 INC lNPUT=2Qi, T APE5=INPUT, OUTPUT* T APE 6 = OUTPUT) 

COMMON /P ARAM/ PHI, PH 10 , OPHI , P3F0 ,PbVS , PND, GEE ,G , NI « NIPi , HH , 

1 NQ,FQ<13) ,NELEV , XELE V < 1 0 1 , Z ELEV ( 1 0 ) ,ELEVCi0l , 

2 ELEV3 C 1 0 1 ,D,H,FNU3,ft W , AQ , XTOT ,OX , DZ, EJDTti, Q0O T , 

3 DIAF, A(7 f laOOIfXUdfil tZCMIOl ,PQS,PHIF,PF, 

4 NC AS£ , L AS TPHI , SA VPHI C 1 C 0 0 » , I SLF PR M , I F AI L 
COMMON /C PROP 3/ XMW , SHRtf , HF G ,PV, RHO L , RHOV , VISL , V IStf , XKL * ST , TF 

OO 500 K=i,10 
CALL rNPT 
11 = 1 
QOT 3=1 • 

000ta=,5 

CALL INGRTNCP3FI 
00 404 1=1,40 

IF( IFAIL • EQ . 1) GO TO 303 
I F( p.jr-PBFG » lu 1, 2C2, 301 


U G U j ? '3 

C C C G 31 
G GO j G 5 
G CO C 37 
G C 3 0 41 
C DO 045 

0 00046 
0 09 052 
0 00 054 
0 QO 056 
0 00 060 
0 00 064 
C GO 066 
C 00 072 

uuO J72 
0 uG 075 
0 CO 10G 


ICi IF* II. t.Q • 1 ) lUJT£> = 2. *0QiJT3 
IFUI.EQ.2) DCOTl;=. 5*09 0 TfJ 
Q D T b = Q 0 T 4 * Q Q D T 9 
• Caul I ;NGR TN (P3F ) 

IFCP3F.GT.P9Fj) 11=2 
GO TO 395 
C 

3D 3 IF ( 1 1 * EQ • 2 1 000 Td= . 5*0Q0T9 
aOTO=QOTd-OQDT3 
IF3Atf=rFAlL 
CALL INGRTNCPdF) 

IF( I^SAV-IPAIL .£1. 1) 11=2 
IFCIFSAV ,E9. 1) GO TO 305 
IFCP3F.LT.PBFu) II = ? 

3G5 CONTINUE 
C 

395 Tl$T~ H JS (PQF-P3F0) /P9F0 

IF ( TEST . LT • O.Oul) Go TO 2C2 
IFCDQOTB .LT. 0.G0i*QDT3l GO TO 202 


4^ 

ro 



RUN X COMPILER <VE*.2. 5M> 


03/23/74. 15. 31.43 


«1ATN 


U 03104 
0 C 3 10 5 
0 00106 
0 00 107’ 
000 110 
0 00112 
oca 112 

0 CO 114 


404 CONTINUE 
202 CONTINUE 
CALL WRT 

IFTNCASE.EQ.OI GO TO 501 

500 CONTINUE 

501 CONTINUE 
ST OH 

ENO 



RU NX COMPILER ( V£<.2. 3M> 


03/2 3/74. 15.il. 43 


00000^ 
CG00C5 
0 1 0 C05 

oooaos 
C00CC5 
0 C0G05 


C CO 0 35 


OOCOC5 
0 00007 
0G0C11 
D 00 013 
0C0C15 
0 00016 
0 00 J31 
0 00034 
0 00 0 41 
0 00041 

C 00055 
0 00C57 
0 00066 
000074 
0 00077 
0 00102 
000105 
0 CO 110 
000113 
COO 116 
COO 125 
G CO 126 
0 00132 
C CO 134 
C GO 136 
G GO 150 


SUBROUTINE OERlV(yP,Y) 

COMMON /FI RAM/ PHI ,PHIQ , DPMI , PQF C ,P 8 VS , PNO , GEE, G, NI , NIP1 .HW, 

1 Ntf *FQC101 ,N£LEV,KELEV(10> ,2ELEV<1Q> tELE V <101 , 

2 ELEV3U0I ,D,H#FNJB,AW,AB,XTOT, OX,OZ,QOT3,QOOT, 

3 DIAF,.A<7,iOQD 1 .XQtlCI * ZQ ( 10 ) * P3S , PHI F , PF , 

4 NCASE,LA3TPHI,SAtfPHI (1C00I * ISL^PRM* I^AIL 
COMMON / CPROP 5/ XMW, SHRV * HFG *P¥,RHOL »P.HOV, VISL, VIS V, *K<- ♦ 31, TF 

t 

DIMENSION YPC41 , YC4) 

YP 1 1 > = 1 . 
po=y<?) J 
N = NEL£V'-i 
00 101 1=1, N 
z=r<i» 

rFcz.c-e*z£LEvf n .and. z.le.zelevu*ih tmk=i 

101 CONTINUE 

DELE V= lELcV-H INK + 11 -tL£V3 Cl MK1 ) / <Z£LEV < 1 MK + i ) -ZEE EV CIMKl ) 

H:l = ELEVi3 (IMO ♦ ( £ LEV3 C IMK + i ) -ELEVBCIMKl ) MZ-7ELE V ( INK) ) / 

1 (ZEbEVCIMK+ll -ZELEVCTMK) ) + 1, 

A 

PHI SPG = PHI Q ’ 

IFCPHIO .EQ # 1,0) PHISP0=C P BS/H)/(1#+PBS/H) 

PH I $ P= (P IS/ H ) /(H3+P05/H) " 

I c ( G £ F .lD. C.l PHISP = 1. 

IF< I3LFPRM .£0. 01‘PHISP*1. 

PHI = Pd S/ (PB+P3S) • 

IFCPHI .LT. 0.11 PH 1=0.1 

IF ( PHI .GT. PHIS P> PHI =PHI S P 

IFCPHl.GT .PHIS J 01 PHI=PHISPG 

IFCPHIF , £Q . 0. .AND. PF .£Q. 0.) PHl=PHin 

IFILASTPHI .£3. C> GO TO 105 

IZ=INTC(ZF1.C1»OZI/OZ» 

PHI=SA\/PHI(H1P1 ) '■ 

IF ( IZ .EC. N IP 1 I GO TO 11$ 

PHl = SAVFHmZ»<-iSAVPHKIZH> *3/ VPH1J IZ) 1 M 7- t IZ - 1) ♦ DZ) /U7 
105 CONTINUE 



RUNX COMPILER IVE?.2.3M) 


08/23/74. 15.31.43 


QERIV 


0 CO 153 
0 00 153 
0 uQ 15Q 

003 205 
000213 
Q GO 214 
0 00 216 
0 00 221 
0 00 223 
0 00 226 
0 00 231 
0 00 235 
0 00 237 
0 00 246 
000247 


PERM=( 3./8.»*Q*< PHI/ ( 1-PHI > ) / 

1 (4. *PHI/ (4. * (l.-PHI) -(l.-PHn**2-2.*ALOG(i.-FHT>-3. ) 

2 3./ (ALOG(l.-PHI) Mi#-(i.-PHI)**2J / (1.M1.-PHI1 * *2 1 » ) 

YP(2l=FNU9*tr ( 31 ✓ iPE RM* 43 > tOELE V*C-EE 
DI5T-0. 

00 112 1=1, MQ 

aisT=oisr*zam 

IM< = I 

IF(Z.LT.OIST> GO TO 113 

112 COMTINUE 

113 YP( 3l=-QDI3 4 Fi( IMK> /ZQ(IMK) 

YP(4)=PHI 

IF ( P3 .C-T. l«00ai*PJS*(l«-PHII/PMU I c A I L= 1 

RETURN 

EM 3 


4k 

cn 



RUNX COMPILER (VER. 2 . 3 M) 


03 / 23 / 74 . 15 * 31.43 


0 LO 003 
0 CO 0 03 
Q uO U03 
0C0003 
0 00 003 
000003 
0 00 0 03 
0 00 003 
GOCOQ4 
000006 
0 coon 
0 00 012 
0 00 013 
0 C0C15 
C03C24 
G0003C 
0OC0 32 
000034 
0 00 036 
0 00 040 
0C0C41 
Q OQ 044 
0 00 046 
0 00 050 
000C54 
0 00 056 
0 00 062 
0 CO 065 
0 uO 06 7 
0 CO 372 
tiOOH75 
0 00076 


SU3R0U T I NE IN6RTN{P3F> 

COMMON /PARAM/ PHI , PHI 3 , QPH I *PBFQ * P B VS* PNO , GEE , G , Nl * NI P 1 , M W * 

1 NQ* FQ ( 1 0 > ,NELEV*XFLEV (l*i> ,ZLLtV (10) , ELEV (101 , 

2 ELHVT ( 101 ,a,H,FNJ9,AW,AB*XT0T,0X,nZ, DOTS* JQOT, 

3 OIAF, A (7, lOOu) , XQ(1G) ,ZQC10) , PBS , PHI F , PF , 

4 NCASE,LASfPHI,SAVFHI C1G00) * I$L C FPM* I FAIL 
DIMENSION Y(4>, YP<4> 

COMMON /CPROPS/ XiA , SHRV ,HFG * P V * RHG L , RHO V. VlSL , VISV, XKL * 51 , TF 

A< 1*1) =0. 

A<2,1)=1. 

IF ( GEE • EQ • 0.) A<2,1I.= P9VS 
A ( 5 , 1 ) - 0 * 

A(4,1)=Q. 

A (5 * 1 ) =PHIC 

l * (PHIQ .EQ. 0.) A (5, 1) =(P J5/H) / (1. +P3S/HI 
IF1LASTPHI ,L 1. 1) A(5*l) -SAVPHKl) 

Y(t)=A(i*l) 

Y (21 =A (2,1) 

Y ( 3 I = A ( 3 , 1 > 

Y ( 4 ) ~ A ( 4 » 1 ) 

IFAIL-0 

CALL DERI V ( YP * Y ) 

NIP i=Nl+ 1 
DO 120 I=2*NIP1 
CALL RUNGS ( Y * YP * CZ) 

DO 11C J= 1*4 
AM t D = Y(J) 

A ( 5 * I ) =PH I 
110 CONTINUE 
120 CONTINUE 

P8F=A(2,NIP1) 

RETURN 

END 



RUNX COMPILER (VER.2*3M) 


36 / 23 / 74 , 1 5 . 31 . 42 . 


SU3ROUT1KE INPT 

CO-1 MON /F ARAM/ PHI , PHI -J , OPHI , Pli^O ,Pfl VS ,PN3 , GHE , G, NT , NIP t , HW , 
x NO,F-ma» f NZLttf,xcLr;wa6)*2CLttf cic»,e:ld/uoi * 

2 EL£\/t)(iQI ,0,H ,FN'JE,A W, AB,XTOT,OX, DZ, SOT T, OJOT, 

3 ' DIA C , A( 7,1 OGui, X H1G) ,ZU(iX ) ,PHS, D H IF, PF, 

4 NCASE,LASTP4I , SJWPHT UGOQ) , ISLFPRM, I FAIL 
COMMON /C PROPS/ X 1 ini , SHRV , HFG , P V , PHOL, KHOtf, VISL , V/ IS*/ , XKL , ST , TF 
DIMENSION 431(61,492(6) 

i} A r A H9l/6*1H / * H02/ 6*1H / 

DATA H/l « T6 5/ , 0/ 2 • 6 3/, R / #60 / , G / 9 * ^ C / 

C 

NAMt LI ST /G 40 AT A/ LIO, TKFLtf N,RHOL,S T,VISL, HFG,6EE, A W,HW, 3IAF, 

1 PH 10, PH IF, OIAVS,UOOT,Na,FQ,NFLr\/,XELEV ,tL fc V 

2 OX, HDl,H02 f X0 r NCASE f LA5TPHI, PF, I3LFPPM 
u 

REA CH 5 * GP3 AT A) 

TR4NK=i. ea*TK£LVN 

CALL PRuF i< HQ, TPANO 

IFUIQ.EQ.OI GO TO U 

ST-4.44H*3.25i*ST 

RHOL=. 4 5 36 *3. 2 U**3 *RHOL 

VISL = 4.44tfM. 2-41**2*36 0 0. *7 I SL 


000023 

HFG= ( 1 35 5 • *2 

• 205) 4 HFG 

0 C 0 026 

13 CONTINUE 


U CO 026 

WRI T E ( 6 , 9 0 Q ) 

HJ1 ,HU? 

GC0036 

WRI T E ( 6 , 9.C l ) 


0 00 3 42 

WRITE(6,932) 

LIQ 

3 00 050 , 

WRITE(6,934) 

TKEL7N 

0 00 056 

WRITE ( 6, 9361 

RHOL 

0 C J 364 

WRITE (6,9031 

ST 

0 00 372 

WRITE (6,9101 

VI S L 

C 00 10 u 

WRITE (6, 912 > 

HFG 

0 CO 106 

WRITE (6, 9141 

GEE 

0 CO 114 

WRITE (6,916) 

AW 

G CO 122 

WRI TE ( 6, 913) 

HW 

0 CO 130 

WRITE ( 0,9231 

31 AF 

& 00X36 

Vi R I T E < 5,922) 

PHI r 


0 00C02 
300302 
0 0D0o2 
0 30 0 02 
0 00 002 
C 00 0 32 
003002 
000002 
000002 

u 00002 

Cuoao2 
0G0 002 

0 GO 002 
C Cu 005 
C COCO? 
0 uOGll 
0 G0C12 
0 00 G 14 
0 uG 017 



RUNX COMPILER. ( VHR. 2. 3M> 


INPT 


03/23/7s. 15.31.43. 


C 00144 
GCatF2 
C CO 16G 
0 CO 166 
C CO 174 
COO 202 
0 03204 
000216 
0 00 221 
0 03 227 
G 00 231 
Q CO 243 
0 CO 246 
COO 254 
0 00 262 
0 GO 273 
0 00 276 
003 276 
0 CO 276 
t C J 276 
000276 
0 00 276 
0 GO 276 
0 C O 276 
0 GO 276 
0 00276 
0 00276 
0 CO 276 
G GQ 27b 
0 00 276 
0 GO 276 
G PC 276 
000 276 
0 C0276 
0 C0276 
C OC 276 
0 GO 27 6 
Q C0276 


5 



DCG 


WRITE(6,924> 
WRITE <6, 925) 
WRITET6,9?6> 
WRT TF (6,928) 
WRITE (6,930 I 

00 6 1 = 1, m 

WRITE (6,932) 


PH IF 
PE 

JI AVS 
Q90T 

m 

I, XG(I ) ,Fiin 


CONTINUE 

WRITE (6, 9341 Nt. LEV 
□0 6 1=1 , NELEV 

WRI T£ < 6 , 93 6 ) I,XEL£V(I> , ELF V ( l ) 

OCR TI N LIE 


WRITE (6,938) OX 

WRITE ( 6, 940 ) NGASE 

WRI T E < 6, 942 I L AS T p H I 

WRI TE ( 6, 944 ) ISLFPRM 

FORMAT (1H1,4X,6A10,/5X,6A1J > 


901 FORMAT {//5X, 

1 5 X 

902 c OR MAT ( 13X, 
904 FORMAT ( 1 jX , 

1 1 12 • 
90 6 FORMAT ( 10 X, 
1 £12. 
JC8 FORMAT ( luX, 
1 E 12 * 

910 FORMATi 10X, 

1 E12. 

912 FORMAT ( 1GX, 
1 L 1 2 • 

914 FORMAT (/1QX, 

1 E 12 • 

916 FORMAT (/10X, 


33H INPUT VARIABLES AND FLUID PROPERT IES S// ) 

, 3 7 H INPUT VARIABLES AND FLUID PROPERTIES*//) 

42HLT0UI0 NJM3ER LIO = 

4 2H T £ MP E R A T U R£ .... T KELVN = 

5 , 16H DECREES KELVIN) 

42HLI0UI0 DENSITY • • * . . RHDL = 

5 , 1 6 H KG/C'J. M ) 

42H SURFACE TENSION ST = 

5 , 1 6 H NEWTONS /M > 

42HL I 3UI 0 VISCOSITY VTSL = 

5,1 NEWTON-SEC/SQ. M ) 

42HLATENT HEAT*,. HEG = 

5 , 1 6H JOULES/KG ) 

42HC-RAVI TATIONAL ACCELERATION.... GE r = 

5 , 2 0 H STANDARD GRAVITIES) 

42HWICK AREA AW = 


1 £12.5, 16H SO. CM ) 

91* FORMAT ( 12X,42HWIC< HEIGHT HW = 

1 t 12 • 5,16H CM ) 

i2C FORMAT ( 1 0 X , 4 2H W I CK P I6ER 0 I AM t TER OIA^ = 

1 . £ 12 . 5 , 1 6H CM ) 


12 ) 



RUNX COMPILER <V£R.?.3M) Q:i Z23/74. 15.3i.V3. INPT 

0 CO 276 92? FORMAT ( lOX , 42H INITIAL POROSITY *HI0 = . 

G0Q276 1 E12.5I 

SCO 276 924 ^CRMAT ( 1UX , 42H FINAL POROSITY..... PHIF = , 

0 CO 276 1 F12.5* 

COO 276 925 FORMAT </ 13 X , 42H F INAL STRESS PF = , 

0 GO 276 t 112.5,166 CM LTQ. ) 

QC0276 926 FORMAT (/ 1 0 X » 42H V A^QR-SP ACE UIAMETER ......... . OlAVS = . 

0 C 0 276 1 El2.5tl6H CM > 

000276 92d FORMAT (/ 10 X , 42H NCMI N AL HE A f - TRANSFF R RATE.... JOOT - , 

000276 1 F12.5,16H WATTS > 

C 00 276 930 F ORMA T ( ljX,42HN0. HEAT-INPUT SECTIONS. NO = , 

C CO 276 1 I?) 

0 00276 932 FORMAT ( 10X, 16MSECTIGN NU M3 E R ,12/ 

0 CO 276 1 15X, 3 7H SECT ION LENGTH <U. = » 

0 00276 2 1 12. 5,1 Q F CM / 

0 00 276 3 15X, 379 HEAT-IN PUT FRACTION FO = , 

000276 4 112.5) 

0 C 3 276 934 FORMAT (/ ICX , 42HNQ, ELEVATION PGINTS NELFV = , 

0 00276 1 T2> 

000276 936 FORMAT (1CX, 21HELEVATION POINT NO. ,12/ 

C 00 276 1 15X, 37HOISTANCE TO POINT XELFV = , 

000276 2 E12.5,10H CM / 

0 00276 3 15X , 37HELEVATION OF POINT F.LEV - , 

0 CO 276 4 L 1 2 • 5 , 1 0 H CM ) 

0 00 276 93 6 FORM AT ( / 1 j X , 4 2 ! t INTEGRATION STEP SI/ E CX - , 

0 00 276 1 €t2.5 t 16HCM I 

0C0276 940 r ORMAT (✓ 13X , 4 2H ANOTHER CASE lO^NO* 1=YLS).*.. NCASE = , 

000276 1 12) 

000276 942 FORMAT I 13X,42HUSE LAST POROSITY 0I5TN LASTPHI - , 

000276 1 12) 

0 CO 276 944 FORMAT < 10X ♦42HKH3JIRE SuLF PRIME-. I3LF PC \M = , 

u LG 276 1 12) 

Q 

0 00 276 PN:J-RH0L*G*(HW/1QG. ) 

0 DO 302 IFTOIAF .£Q. 3.01 GO TO II 

0 CO 303 P9S=3.?46*R*H*$r/<<;UAF/103.)*PNn> 

0 CO 311 II r p (G£E.NE.O.) GO TO 12 

0 CO 312 t>3VS-<4.*3T / CDIAVS/100. n/PN9 


■p. 

V£> 



RUfMX COMPILER IVER.2.3MI 05/2.3/74, 15.31.43. 


0C0 317 


12 

CONTINUE 

0 CO 317 



XTOT-O. 

0 CO 320 



90 16 1=1, NQ 

0 CO 322 



XT0T = XT0T4-XQ(I) 

0 CO 325 


16 

CONTINUE 

0 00 330 



A9=AW/HW**2 

Qufl 333 



DO 13 1 = 1, NEL EV 

0 00 335 



ZEL2V(II=XELEV(II /XTOT 

0 00 343 



ELEV8 ( II =ELE V ( II /HH 

0 CO 343 


18 

CONTINUE 

0 00 346 



00 23 1=1, NQ 

0 00 350 



ZQ(I)=XQ(II /XTOT 

0 03353 


20 

CONTINUE 

0 GO 356 



NI=XTOT/DX 

0L0 361 



0Z=0X7XT0T 

u GO 363 



IF ( PHI 0 . ME . 0.01 GO TO 50 

0 00 364 



IF1PHIF • ME • O.G) GO TO 50 

0 00 365 



IF ( PF .NE . 0 .01 GO TO 50 


u 

r' 

c 

_ 

OP T I MUM-HOMOGENE QUS- WICK CALCULATIONS 

0 G 0 366 

V 


IF (01 AF • N£ * O.QI GO TO 30 

0 DO 367 



PHIG = 0.79 

0 C O 371 



0 BS = 2.MELEVd(NELEV)+i.)*PHIG/(l.-PHIC) 

0 GO 400 



01 AF=1 00 . *3 . 246*R*H*ST/ (P3S*PND> 

C Q 3 405 



WRITE (6, 950 I OIAF,PHIO 

C CO 415 



GO TO 50 

0 00416 


3 0 

PHI Q =P3S / ( 2 * MEL E VO < NELE VI *1 • ) +PBSI 

0 lQ 423 



WRITE (6, 9521 PH I 0 

C GO 431 


95 0 

FOR HAT (///10X»4QH OPTIMUM HOMOGENEOUS WICK* FIOER 

0 00 431 


1 

E 1 2 • 5 , 16H CM, POROSITY =,£12.5) 

COG 431 


952 

FORMAT (///1GX,47H OPTI MUM- POROSITY HOMOGENEOUS WICK 

0 00 431 


1 

£12. 5) 

0 00 431 


5 0 

CONTI NUF 

0 GQ 431 



FN Jb=( VISL/RHOlI MQDOT/HFSI* (XTOT/iC.O. > / 

£ GO 431 


1 

(PNOM 014*/ 130. >**2MHW/1Q0.I**2) 


o 


IN PT 


OTA. ' = , 
i °0R05ITY 


ur 

o 



RUNX COMPILER Ci/a«2.i M) 


0 i/23/74. 15. 31.43 


• INPT 


0 00 446 


C CO 460 


0 CO 473 


0 GO 501 

COO 514 
0 00514 
000515 


C CALCULATION OF FINAL STRESS PSF D 

C 

C 1. PF SPECIFIED IN GRADED- P3R0SIT i DESIGN S 

C 

IFtLASTPHI.EO.O .AND. PF.GT.0,0) pbfo=pf/hw 
C " 

C ?. PH IF S P l C I F I ~ i) IN GRADEO-POROSI TY DESIGN* 

IF( LASTPHI . EO.O .AND. PHI^.GT.O.O) PB^C “P3 S* < 1 • -FHI F) /PH IF 
C 

C 3. *BF0 SET BY PHIF OF LAST CASE: 

n 

sJ 

IF (LASTPHI .EQ.l ) P 3 F Q =P 3S* ( 1 . - S AVP H II NIP1 )) /S AVPHX C NIP1 ) 

C 

C 4. HOMOGENEOUS WICK - P3F3 SET BY PHI3 i 

r* 

U 

IF (LASTPHI. EO. 3 • AND • PH IF+PF . HQ . C J P3F0 =PBS* (i.-PHIi;)/PHI3 

j "\ 

o 

2 2 CON 1 TN UE 
RETURN 
END 



RUN* COMPILER <\1 tR.2.3M> 


03/23/74. 15.31 .43. 


C 

c 

/> 

0 OQ 0u5 
Q GO 005 
0 C 0 u 0 5 
000005 
000005 
0 CO Q 05 
000005 
0 GO 005 
0 00 005 
000005 
OCQQ05 


SUBROUTINE PROPS <L» Tl 

THIS ROUTINE COMPUTES FEUIO PROPERTIES FROM DATA FITS 

COMMON / CP R OPS/ XMW ,-SHRV ,HFG ,PV ,RHOL ,RHOV, V'I SL , VISV, XXL * ST* TF 
31 MENS ION A11(7I , A2K7I * 

1 A 31 (7), A32<7>, A3 5(7), A34(7), A35(7), 

Z A 4 i ( 7 1 t A42C71* A4 3<7>, A44C7), A45(7), 

3 A 5 1 C 7 A52(7), 453(7), A54(7I, A55(7>, 

4 Ab 1 ( 7 ) , 462(7), A 6 3 ( 7 ) , A64(7>, A65C7), 

5 A 7 1 ( 7 ) , A 7 2 ( 7 ) , A73f7l, A 74 ( 7 > , A75<7), 

6 A 3 1 ( 7 ) , 40 2 ( 7 1 , 433(7) , A J4 ( 7 ) , 4 35(71* 

7 A 9 i ( 7 ) , A 9 2 ( 7 ) , A93(7>, A 94 < 7 ) , A 95 (7), 

3 A 1C 1 ( 7 ) , A 10 2 ( 7 ) * A103(7|, 4104(7), Ai05(7), 

9 A11K 7) , 4112(7), All3(7) t AH4(7), A115(7) 


C WATER (32F<T<400F) 


CC0G05 

ja r a 

A 1 1 < 1 1 , 

DU0 005 

¥ 

491.7, 

0 00 GG5 

DATA 

A o 1 ( 1 ) , 

0 GO 005 

♦ ' 1. 

3555636 , 

0 00 005 

DATA 

A4 1 ( 1 ) , 

fiG00G5 

♦ 1219. 5 508, 

0 00 0 95 

DATA 

A51 < 1 > , 

D 00 0 35 

* 14 

.199322, 

0 0Q 005 

OAT A 

A6 1 ( 1 ) , 

0 GO 005 

♦ 58 

•481 766, 

QC0 105 

DATA 

A 7 1 (1 ) , 

C GO 005 

* 7, 

4432132, 

0 GO 1C 5 

DATA 

ASK 1) , 

0CQ005 

* 52 

.625735, 

0 00 J 05 

DATA 

491(1) , 

0 CO G05 

* -11.6a486, 

0 GO 0 0 5 

DATA . 

A1C 1 (l) , 

0CQ0C5 

* -1. 

0535655, 

0 CO GO 5 

DATA 

4111(11, 


A21 ( 1) / 

13.016/ 

A3 2 ( 1 > , A i 3 ( 1 ) | 

4.9575765-5, D. , 

A4 2 ( 1 ) , A 4 3 (II , 

5* 70551 5E-2,-4. 454538 E “4, 

A 5 2 ( 1 1 , A53 1 1 ) , 

-6.5267262, -.31013069, 

Ab2(l), A6 3 ( 1 ) , 

2. 586296E-2,-3. 54721 2 £>5, 

A7 2 ( 1 ) , V73(l> , 

-6.Q175647, 7904 0982, 

482( 11 , A33 { 1 ) , 

-.26276095, 5. 0 3327CE-4 
A92 ( 1 ) , A 9 3 ( 1 ) * 

i* 1341087, G. f 

A 1 G 2 U ) , A 1 j 3 ( 1 ) , 

, 352 2 992E- 3, -6.44 605 02-6, 2. 
A 112(1) , 4113(11 , 


A 34 ( 1 ) , 

A 3 5 ( 1 ) / 

0 ♦ , 

0./ 

A44 { 1 ) , 

A45 (1) / 


0 ./ 

A54 ( 1 ) 

A55 (11/ 

0 . , 

1 ./ 

A 6 4 ( 1 ) , 

A65 ( 1 ) / 

u • , 

3 ./ 

474(11 , 

A 7 5 ( 1 ) / 

1 . , 

0./ 

A 3 4 (1 ) , 

A 8 5 ( 1 ) / 

4 11 82 35- 7 , 

1. 469625-10/ 

A 94 ( 1 ) , 

A95 (11/ 

- ■ ? 

0./ 

a i a 4 ( 1 1 , 

A1C 5 (11/ 

15232 7E-4, 

ii . / 

All 4( 1 ) , 

4115(11/ 


en 

ro 



RUNX COMPILER (VER.Z* 3M) 


33/23/74. 15. 31.4 3 


PROPS 


C GCC05 


G 000G5 
D 0 G CDS 
G 00G95 
L CQ 005 
G GOG05 
CC0GC5 
G GO CO 5 
0 G G 0 0 5 
C Q 0 ii 0 5 
GOO 005 
DG0005 
0 CO 0 35 
0 CQC05 
G GO 0 05 
00^005 
G u O 0 05 
C GO uC5 
G GO 0 35 
0 GD 005 
C GQCC5 


0 GO 005 
0 00 005 
0 09 305 
0 0GC95 
0 00005 
0 00 005 
0 GO GC5 
0 G 0 u 0 5 
G CG005 
C GOOC5 
G CJ005 


♦ - i.43 7350E-3, 9 . 7 1 722 3E-5 , -2 • 2 3T75 7E- 7 , 2. i 171 95b- 1 u 


AMMONIA ( - 1 U 7 . 9 F < T < 1 9 0 F ) 


DATA 

All( 2) , 

421(2)/ 


351*3, 

17*032/ 

DATA 

A 3 1 < 2 ) * 

A 3 2 ( 2 ) * 


1.31, 

3 • i 

DATA 

A41(2# , 

A4 2 1 2 > * 


* l*Q93251E+3»-2*4b2955E*0, 

DATA A51(2) ♦ A52(2> , 

* 1. 39 23 74 E* 1, -4. 921 740 E*0, 

DATA A 61 ( 2) , A62(2I , 

* 7*343766c. + l,-l* 17240 54-1, 

DATA A 71 { 2 > , 472(2) , 

* 1.266966E+1»-1*113379E+1, 

DATA A81I2) , A32<2> , 

* .*» * 5 3 7 C 46E ♦ 1 * - 2 • 496 424E- 1 * 

DATA A 9 i ( 2 ) * A92 12) » 

* - 7 • 0 7 3 3 9 6 E + 3 * 1 • 9664 9 4E+ 3 , 

DATA A 1 G 1 < 2 ) * A132(2J, 

* -4*16 01864-1* 3 #944713 E - 3 , 

OAT A A11K2) , A 112 ( 21 t 

* 6*4265 DIE -3, -7*6046 4 1F-6, 


A33 <21 * A34 (2) , 

C • * 3 « ♦ 

A4 3 ( 2 ) * A44 ( 2 ) * 

4. 97643 OE-3 *-4.4749675-6, 
453(2) * A54 ( 2 ) * 

2*0 65018 E-lt-7 #5795 9 7c-2* 
A63(2>, A64(2), 

1,3 3170 7c- 4 i-l *649 9136- 7, 
A 7 3 < 2 ) , A 7 4 ( 2 ) * 

2. 99312 b£*u» -4* 639769E-1* 
Ad3<2) i A 64 ( 2 1 , 

6* 6231664-4* — 7 • 941 805 r -7 * 
A9 3 ( 2 ) * A 9 4 ( 2 ) * 

-4.728715E+2 f 5.G54C66E*lt 
A1G3<2>* A 1 G 4 1 2 ) * 

-6.537242E-6* 3 . 0694354-9, 
All 3 ( 2 ) * A114<2), 

-7.6997595-9* 3. C23573E-12, 


METHYL ALCOHOL <-140F<T<38GF) 


DATA All Ui , A21 < 3) / 

* 322*7* 12.042/ 

DATA A 3 1 ( 3 ) « A32(3I * 

* 1.20 3, 0* , 

DATA A41< 3) * A42 1 3 ) , 

* 3* 790 5 4bE *2,-2 . 4781Q 5£ + 3 , 

DATA A 5 1 ( 3 ) * A52U), 

* 1.5c 541 1E + 1 1 - 9. 24 0 63 0£ + J * 

DATA A61I 3) * A6213J , 

* 1. K7d :, 3e>i, 2 . 9 3 2 6 3 6 E - 1 1 

DATA A7.1 < 3 ) , A72( 7) , 


A3 3 (3 ) , A 3 4 ( i ) , 

C * * 0 * * 

A43 (3) , A 4 4 (3 I , 

6 • 4 1 & 62 9 £- 3 * -7 * 0 u 4 1 9 5 E - 6 * 
A53(3) , A 54 13) » 

3 . 3 661 7 6£> G , - 1 . 969 2 0 0£ + G * 
A6 3 ( 3) , A64 ( i) , 

-3, 417672 E -4, 9*7612629-7, 
473(3) * A 7 4 t 3 ) * 


7 * 53 0 81F - 1 4/ 


A35 (2)/ 

. 0 ./ 
A4 5 < 2 ) / 
0 • / 
A55 ( 2>/ 
0./ 
A&5 ( 2 ) / 
0./ 
A 7 5 1 2 ) / 
0 ./ 
A85 (2)/ 
.552154E-10/ 
A95 ( 2 ) / 
2. 024369 E + n/ 
Ain5(2l./ 
0 ./ 
A115 (2)/ 
n . / 


A 7 5 (3)/ 
n ./ 
A45 < 3 ) / 

2. 2144393-9/ 

A55 (7)/ 

3. 7 6 965 6 c- 1/ 

A65 (7)/ 

4. TG 5 32E - 1 n / 

A 7 5 ( 3 ) / 


cn 

Ca> 



RUNX COMPILER (VE0.2.3M) 


Di/23/74. 15.31.43 


PROFS 


GOO 005 
C'CQ0n5 
0 03 C05 
3 03 0 05 
0 CO 005 
003005 
0 CO 005 
0 00 305 
GOG 005 


0 03005 
0 03 005 
000005 
000005 
0 CO 0 3 5 
0 30 005 
d CCCG5 
0 G 0 0 0 5 
0 C03C5 
C 0QQC5 
0 GQ305 
CGQG05 
0 Cu 0 0 5 
0 CO d ti 5 

C GO 0 05 
b CO CC5 
Q 3 u 0 05 
C GO 305 
000095 
U OG 0u5 


C 03 035 
0 CO 0C5 
0 to Q05 


* 1. 59 U64E+1 *-2 .10 9398tc + l, 1 . 1 4432 6E+ 1 » -4 . 2 786 4 3 1> 0 . 

QmTA Aftt(3), 482(3), A83C3I* 184(3), 

* 2 • 2 8 5 3 63E + 1 , - 1 * 1 5 3 1 6 9f. - 1 * 2. 30 3795E-4 , -2. 1551 27£- 7, 

QAT A A 9i ( 3 ) , A 92 l 31 , A9 3 (3 ) * 194(3 1 , 

* l»522596E+2,-l»?38769E+2» 3. 113 34 4L+1,- 3.31341 Oic+O, 

SATA A 1 0 1 ( 3 I , A 1 0 2 ( 3 1 i AID 3 (3), A1C4(3>, 

* 3.944433E-2, 2 . G97417E-4* -b. 71 C 59 7n- 7 , 7. 394364E-10 , 

SATA A111C 31* A 11 2(3>* A 11 3 (31* A11^<3I. 

* 5*793525 £- 3, -1.4u4494E -5, 1*205 62 0E> 8*3. 5 166 2 9 c *'12* 


FREON- 21 ( -55P<T<305 C I 


DATA 

All ( 4) , 

A 2 1 ( 4) / 

* 

248.7, 

132.93/ 

DATA 

A 3 1 ( 4 ) , 

A32 ( 4) * 

* 

1.175, 

A. , 

DATA 

A4l(4) , 

A4?(4) , 

* ft. 

537625E+1, 

4. 6 3 655 8 £-1* 

OATA 

A51 ( 4) , 

A5 2 ( 4) , 

* i. 

67.732E+0, 

1.573170E+1* 

DATA 

A61(U , 

A6 2 ( 4 > , 

+ 1. 

332756E+2 , 

-3 * 261 /57E-1 * 

OATA 

A 7 1 ( 4 ) , 

A72 ( 4) , 

♦ ft. 

5 3 4 322E + i * 

-1.6625 7 5£ + 2* 

OATA 

A 3 1 ( 4 ) * 

A ft 2 ( 4 ) , 

* -ft. 

34747^6+ 0, 

ft • 5 3*1 1 1 oc* 2* 

OATA 

A91 ( 4 ) * 

A92 ( 4) , 

* -1. 

83353 3 E+3* 

1.199366E+3* 

DATA 

A1C 1 ( 4) * 

A iu 2 < 4 1 , 

♦ 

75019 9 E-l* 

- 2 . j 3 5 4 ft E - 3 * 

DAT A 

AilK 4) , 

A 1 1 2 ( 4 ) * 

* - 3 • 

2489712- 3, 

4 • 98 4 ft6 9tc -5 » 


ETHANE ( - 1 3 5 F < T < ft 5 F ) 


A3 3 ( 4 ) , A 34 ( 4 ) . 

0 • * O.i 

A43(4>, A 4 4 ( 4 I ■ 

1.6316 8 5 t- 3, 2.C56597E-6 
A5 3 (4) , A 54 ( 4 I 

1. 60 7959541 , 5 .2592436+ Q 
A6 3 ( 4 I , A 64 { 4 ) 

1. 11165 56-3, -1.6 11 7 2 ft- -6 
A 7 3 ( 4 ) * A 7 4 ( 4 ) 

1. 2525? 16 + 2, -4. 265ft ft 26 + 1 
A ft 3 (41 , 464(4) 

— 2 • 7 b 7 3 8 6 t — 4* 3 . 6 4 3 7 2 4 - - 7 
49 3(4* * A 94 (4 ) 

-2.94471 It +2 * 3 • 2 1 5 0 7 6 l-+ 1 
A103C4I , A 1 C' 4 ( 4 ) 

5.713512 t -6, -6. 391 8C 25-9 
A 1 1 3 1 4 ) * A 1 1 4 < 4 I 

-1. 13 374 7 2-7, c. 2356536-11 


OAT A All ( 5) , 
* 161.6* 
OAT A A 3 1 ( 5 ) * 


A21 ( 5) / 

*0.07/ 

A3 2 (5), A3ft(5), A 34 (6) 


5 • ft C 9 9 0 ft E - 1 / 
A 35 (3)7 
7.55 172E- 1 1 / 
A95 (3 ) / 
1, 325C51E-1/ 
A 1C 5 (3)/ 
3.196966-17/ 

a 1 1 5 m / 

ft. 670296-15/ 


A35 (41/ 
C ./ 
A45 (4) / 
-1. 0149485-9/ 
A55 (4) / 
-6. 2095310-1/ 
A65 ( 4 ) / 

f .9 j6674F-ir/ 

A75 (4)/ 
5. I75467t>’/ 
A 85 ( 4 ) / 
-1. 7 6 7 8 7 E - 1 C / 
A95 (4 ) / 

-1 . 315723 c + 0/ 
A 1 06(4)/ 
2 * o5 3 4 OF - 1 2/ 
411514)/ 
-2. 25959E-1 4/ 


A 35 < 5 ) / 


4 * 



RUN X COMPILER 

(VER.2. 3M> 

03/23/74. 15.31 

0 GO 0 05 

* 1.18, 

Or, 

G * , 

C CO 005 

DATA A41 ( 5) , 

A42 (51 , 

A43 <5 ) , 

0 GO 005 

♦ -4.2789342+3, 

4.5732542+1, 

-1.7194812-1, 

0 00005 

OAT A A5 i ( 5 > , 

A52 ( 5) , 

A5 3 (5 ) , 

GOO 005 

* 4.513520E+1, 

-5. 3G 32732+1, 

3.3885052+1, 

0 CO 0 05 

i ) A T A Abl ( 5) , 

A62 ( 5 ) , 

A6 3 ( 5 1 , 

0 00 005 

♦ - J.433C14E + 2. 

3. 3010412+0, 

-l.478827c.-2* 

0 G G 0 9 5 

D AT A A71C5) , 

A72 ( 5 ) , 

A73 (5) , 

O CO 005 

* 9. 3310802*1, 

-1.4637312+2, 

8 . 4 22926 E+ 1 , 

G GO 905 

DATA A 3 1 ( 5 ) , 

A82 (5) , 

A 3 3 < 5 ) , 

0G0CG5 

* -1. 72 39432 + 1, 

1. 931929E-1, 

-7. 953422E-4, 

GC0C05 

data A91( 5) , 

A 9 2 ( 5 ) , 

A9 3 ( 5 1 , 

G 0 0 t G 5 

* 2.9908552+4, 

-2.0174352+4, 

5.035*132+3, 

G GO 09 5 

DATA 4101(5), 

A 1 0 2 t 5 ) , 

Alu 315) , 

0 GO G i 5 

♦ -1.142869E+0, 

1.3170962-2, 

-5.07252 5E-5, 

0 GO G 95 

DATA Alii (5) , 

A 112 ( 5) , 

All 3 ( 5 ) , 

OCOU05 

♦ 1* 1 2 37G9E-2 , 

-8* 3396222-5, 

2.759121L-7, 


C METH4N2 (~ 2 BjF<T<-i 2 nr) 

r 


0 CO 005 

Data ah (5), 

A 2 1 ( b 1 / 

4 

it C 0 i. 9 5 

* 163.2, 

lb. 94/ 


G GO 005 

DATA 431(6) , 

A3 2 ( 6) , 

A3 3 ( 6) 

C GO 005 

* 1.32, 

0., 

u • 

G 00 C 05 

DATA A 41 (6) , 

A42 ( fc) , 

A43 (6) 

G 0.0 005 

* -1.124Q912+3, 

2 . 4 25142E + 1»- 

i. 53930 7E-1 

C GO G05 

DATA 451(61, 

452(6) , 

A53 (6) 

0C0u05 

* 1 * 3656342 + 0 » 

6.55761 f t. + 3 , - 

1. 74.657 OE+C 

C GGC05 

DATA A 6 1 < 6 ) , 

462(61 , 

A6 3 (6) 

000005 

* i. 4699512+1, 

3.8141312-1,- 

3 • 22354 2 E- 3 

0 GO G 05 

DATA A 7 1 ( 6 ) , 

A 7 2 ( b ) , 

A7 3 ( 6 ) 

C 00 005 

* 6.391062E+1,- 

5. 495G63E+1, 

1. 3374932+1 

0 00 005 

DATA A 3 1 (6) , 

A3 2 ( 6) , 

453(6) 

0G9 JC5 

* -3.526486E+G, 

2 • ii 24132E- 1, - 

1. 5408205-3 

GG0G05 

DATA A9116), 

A92(6) , 

A 9 3 ( 6 1 

G GO 305 

* ). 11 36312+3,- 

6* 7o662 9E + 3, 

1 • .8 6 0 5 2 2 F + 3 

003005 

J AT A A 1 o 1 ( 6 ) , 

4102(6) , 

AiC 3(6) 

G 0 G C 05 

* 3* 4864732-1,- 

1 . 72 J 9552-3, 

3,6992972-7 


cn 

cn 


43 


PROPS 


2 . > 0 ./ 

A44 ( 5 I v A45 (51/ 

2* 8404392-4,-1. 7568S9E-7/ 
A54(5I. 455(51/ 

-9*1657732 + 0, 9* 15 4704 2-1/ 
464(5), 465 ( 5 ) / 

2 « 4511 66 2- 5,-1.5186822-8/ 
474(5), A 75 ( 5 ) / 

-2* 1918412+1, 2* 1298032+0/ 
484(5), A 65 ( 5 ) / 

1.35510 3E-6,-6. ; IG50 66>1C/ 
A 9 4 ( 5 I * A 95 < 5 ) / 

-5.697Q99E+2, 2. 5928702+1/ 
A 1 0 4 ( 5 I , A 1 G 5 < 5 > / 

8*3902 942-3,-5* 118602-11/ 
A 114 ( 5 ) , A 115 ( 5) / 

-4* 39 3 4 3E.-1J, 2 *65146 85-13/ 


A34<o), A 35 ( 6) / 

J 0./ 

A44 (6) , A45 (6) / 

4. 547 11C £-4, -4* 90 6714 2-7/ 
A 54 ( 6 I , A 55 ( 6 ) / 

5,3033472-1,-3*2571582-2/ 
A 6 4 ( 6 ) » A65 (6) / 

1* C76 42 02-5,-1* 3531232-8/ 
474(61 , A75 (6)/ 

-2.799C4 8EMJ, 1.56 7595 2- 1 / 
484(61, A 35 (6)/ 

4, 7197152-6,-5* 211675 2- 9/ 
A9 4 { 6 ) , A95 < 6 I / 

-2.3215852+2, 1. 07 440 Or +1/ 
4134(6) i 41 05(6)/ 

1*8407472-8,-3. f 3 32 32-11/ 



R UN X Zi 

dkpiler 

< VE R .2. ?M> 

08/23/74. 15. 31.43.- PROPS 

C l Q 005 


DATA A 1 11 ( 6 ) 9 

A 11 2 < fc) , 

All 3(6), All 4 (6 ) , 

4115(6)/ 

G CO 00 5 

p 

♦ 5.4127976-3,-5 

*4639476-5, 

2. 3562406-7,-7 . 617006-10,8. 

1464956-13/ 


n 

Ly 

r 

NITROGEN (- 

340F< T<-25c 

F) 


oooacs 


0 A r A a 1 1 ( 7 ) « 

A2 1 ( 7 ) / 



C LQ GC5 


•* 113.9, 

28.016/ 



C 090C5 


DATA A 3 1 ( 7 ) , 

A3 2 ( 7) , 

A3 3 ( 7 ) , A 3 a ( 7 ) , 

A35 (7)/ 

D 0 0035 


* 1.40, 

0 . , 

C . , G . , 

0./ 

G0GG05 


OAT A A4 1(7), 

A42C7) , 

A4 3 (7) , A 44 < 7 ) , 

A45 (7)/ 

G GG 035 


* 7.64 89746+1,-2 

.3055566-1, 

5. 3175996-3,-2. 3407156-5, 

C . / 

003005 


DATA A51( n . 

452(71 , 

A5 3 ( 7) , A54 ( 7 ) , 

A55 < 7) / 

C CGGG5 


♦ >. 21 7173c + 1, - 1 

.4 312396 + 1, 

3.0547646+0,-3.1337776-1, 1 

. 17 644 96-2/ 

C 00 0 05 


0 A T A Afcl ( 7) , 

A62 ( 7 ) , 

463(7) , A 64 ( 7 ) , 

A65 (71/ 

G 00 035 


♦ 7.29371 36+1, -3 

. 32 32326-1, 

2.2314696-3 , -7 . 4 73 6 326- 6, 

C./ 

000 005 


DATA A 71 171 , 

472 ( 7) , 

A73 (7) , A74 (7) , 

A75 (7)/ 

0 00 0 35 


* 2. 10 28 C 26 + 1 » - 7 

.5037276+0, 

9.6132736-1,-4. 8611 16 E- 2, 

0./ 

003005 


DAT a A81 ( 7 ) , 

A62( 71 , 

A 3 3 ( 7 ) , a 34 ( 7 ) , 

AA5 (7) / 

0 L G 3 0 5 


* -1.1 30 70 96+1, 3 

.530 9376-1, 

-3. 6809436-3, 1.6373146-5,-2 

» 5C 92 66 c- 3/ 

C 00 0 05 


DATA A 91 ( 7 ) , 

A92(7>, 

A93 (71 , A94( 71 , 

A95 ( 7 ) / 

0 00035 


* 1.71 3 6706+4,-1 

« 3715916 + 4, 

4. 13 39456+3,-5.45 3 61 56+ 2, 2 

. 7166246 + 1/ 

C OOOG5 


DATA 4101(7) , 

A lu 2 < 7) , 

Alu 3(7) , Ai 9 4 ( 7 ) , 

A1 05 (7)/ 

C 0 00 05 


♦ 1.1730006-1,-7 

» 9 924246-5 , 

- 1. 4 3 151 5c- 6, 3 . , 

0 ./ 

010005 


DATA A 1 1 1 ( 7 ) , 

4 1 1 2 ( 7) , 

All 3 17) , fill 4 ( 7 > , 

4115(7)/ 

Li 0 0 305 


* 1.6360316-3,-3 

. 7639396-6, 

-4. 3793716-3, 1. 2703966-10, 

0./ 

C 0GQD5 

u 

c 

IF <L . 60 . 0 > GO TO 20 



0 C u 3 C 6 


T 2 = T*T 




DOG 0 0 7 


Ti = T?*T 




OCOLli 


T 4 - 12*12 




0 C0C13 


TR - 1 3 01:. 0/T 




0 03115 


TR7 = TR + TR 




0 30 317 


TR> = T k 2 * T R 




C CO 321 


TR4 = TR?*TR2 




COG 323 


ALT-ALOC (T) 




G GO 027 


AlT2=AlT»ALT 




0 CG031 


A L T 3 = A L T 2 * A L T 






RU4X COMPILER 


08/23/74. 15.31.43 


PROPS 


C C0 0 33 


ALT4=4LT?*ALT? 


C FLUID PROPERTIES 


000 035 
0 0003/ 
0 00 041 
000053 
t G0C65 
0 00103 
l CO 115 
0 CO 1 33 
£ OC 151 
0 00167 
0 CO 201 
0 00213 
000215 
0 CO 217 


0 L u 2 20 
0 00 220 
C Cu 221 


TF - All (L) 

XMH-A21 <L> 

SHRtf - A3l(U+432lU*TfrA33lL>*T2+A34(l)*T3+A35(U*T4 

HFG = 44i<0+442<L>*T«-A43<L)*T2*A44<L>*T3 + A45<D*T4 

PV = EXP<A51(L> ♦A52(L>*TR+A53iL)*TR2+A54(L>*TR3+A55(L>*TR4) 

PHDL = Abl(LH-A62(L>*T*A63<L>*T2 + A64(U*T3 + A65(U*T4 

RHOV = F.X=>(A7HL) + A72(L)^TRf A73(L)*TR2fA74 (U »TR3* A75<L) *T<?4) 

i/ISL - EXP ( A ail L »«■ *82(L> »T+ A 8 3 (L ) *T 2* A6 4 (L) *T 3* A 35( L > *T4 » 

V/TSV = £XP( A91(i-)+A 92(L) *ALT + A93 (LJ *ALT?fAy4(U *ALT 3 + A95 (L) *ALT4> 
XKL = Aid (U * A102(L> •T*A103 (L) *T2+ a1C4 (L) +T34 A1D5 C'U*T4 
ST = AU1(LH-Ali2a> *T+AU3<L» *T2*A 1 14 (U *T 3+ A 1 15 { L I *T4 
tf'ISW = VI5V/4.1697504F3 
VISL-VISL/4. 16975 3 4E 8 
PF TURN 
C 

2D GOMTINJE 

RETURN ■ " 

en:j 



RUNX COMPILER (VER. 2. 3M) 


03/2 3/74. 15. 31.43. 


0 0QCO6 
G00QC6 
060606 
o oacoo 
£ 63006 
0 uO 006 
0 0C 006 
6 00606 
0 60 0 1«j 

0OGG13 
0C0C22 
0 C 0 026 
0 00 030 
0 00 037 
0 C O 643 
060045 
0 00053 
0 00 057 
0 C0061 
0 60076 
U GQ1C2 
0 60103 


SUBROUTINE RUN5E ( Y, YP, DZ ) 

O OH MON /PAR AM/ PHI, PHI 0, DPH I ,POFG * PB VS, PND , GEE , G , NI t NIPl»HH v 

1 NQ.FliiQ I , NELEV, X£L£V«10 > tZELEy ( 101 , ELEV(lO) . 

2 EL-VBUO) ,0, H,FNUS,AW,AB»XTOT , OX , DO, QOTB , QOOT , 

3 OIAF r AI 7* 100 0 I *XQ(10) , Z Q < 10 ) , P3 S , PHI F, PF , 

4 NCASE^tASTPHI .SAtfPHI (106 0 ) , ISLFPRM.IPAIL 

DIMENSION Y (41 , YP(4» , YP1<4> , YP2C4J , YP3<4», YP4<4I , YOC4I 
COMMON /C PROPS/ XMW,SHR\/,HFG*P\/»KH0L,RHCV*i/I$L,Vl3V,XKL»ST,TF 

0011=1,4 

ypiu>=yp(i> 

1 YO( I) = Y( I) ♦YPK U*OZ/2. 

CALL DtRItf ( YP2, YD) 

DO 2 1=1,4 

2 Y0(I>=Y<I)+rP2<I)*0Z/2. 

CAuL O.IRIVMYP3, YDI 

DO 3 1=1,4 

3 YO( II = Y( I) f YP3( I)*OZ 
CALL D ER I \M YP4, Y C) 

DO 4 1=1,4 

4 Y t II = Y < 1 1 + C YP1< II *2 . * YP2 { I) +2. *YP3< I M- YP4( f) I * ( DZ/6 • I 
CALL CERI\/(YP,Y I 

RETURN 

PNID 


cn 

00 



RUNX COW^IL^R <V£R.2*3M> 


08/23/74. 15,31. 43* 


CCQ302 
0 00 0 32 
0 00 0G2 
0 CO 002 
0 00 0 02 
0 0C u 32 
0 00002 
0 CO C 03 
000014 
0 C0C16 
0 CO 022 
0 00 030 
0 C0C37 
0 00 C 37 
000055 
0 00060 
Q C J 06 3 
0 00 075 
0 uO 103 
0 00 102 
000106 
0 CO 116 
0 CO 123 
0 CO 144 
C CO 147 
0 00147 
0 OC 147 
0 UO 147 
0 10 147 
C CO 147 
CCO 147 
0 03 147 
0 CO 147 
0 CO 147 
D OC 150 


SUBROUTINE WRT 

COMMON /P ARAM/ PHI, PHIQ , DPH I ,PBFO , PB VS, PNQ , GEE , G , NI , NIP1 , HH f 

1 MQ,PaC10) ,NELEV,X£LEV<10I ,2 lLCV(10I , ELF VI 1C) , 

2 tL EV 3 ( 1 0 ) ,0,H,FNUS,AW f A8, XTOTtDX, QZ, QGTtffQDOT, 

3 DIAF, A (7, 13 3C I ,.Xa< 10 I , ZQ (IQ I f PoS , PHIF, PF, 

4 NCASE,LASTPHI,SA¥PHT C1CDC), ISLFPRM, t FAIL 
COMMON / CPROPS/ XMW , 3H Rtf , HFG ,P tf , RHOL * RHOtf * V ISL • V I S'J , XKL * ST , TF 
I M < = 1 

C-M4SS=A < 4 f NIP1I* ( AW/ i 00 0 3 . ) * RHOL* (X TOT/ 1 00 . ) *1 3 0 C • 

00 53 J=1,NIP1 
A(i»J>=XT9T*A(l f U> 

A(2,J) = iCa.*PN ; )*A<2,J)/ ( RHOL *G) 

IF( A ( 1 , j) .G£.XtLE'9( IMKI > IMK=IMK+1 

A 16 , J) =ELEV < I M<- 1 ) I A C 1 , Jl - XELtV ( IMK-ll ) MELEtf 1 1 *K* -E'L£V< IMK-11 ) 
1 / (XELEiHIMKI -XELEV < I MK- 1 ) > +H w 

IF (Gl£ .£■}. 0.0 > A(6,J)=D. 

SAtfPHI(J)=A(5,J» 

A<7, J)=(ll)0.*P.MO/IRrtOL*SI> *P3S*!1.-A«5,J(| /A(5, Jl 

53 CONTINUE 

QOT=QOTd*ODOT 
WRITE <6, 960 ) 

WRI T E (6 , 962 ) OQ T * GM ASS 
00 9 C J= 1 * NIP 1 

WRITE(o, 966) A ( 1 , J J , 41 ?, Jl , A (6, J) , A ( 5, J> , A ( 7. J ) 

90 CONTINUE 

960 r ORMAT C1H1, 4X t 15H c INAL SOLUTION //> 

962 FORMAT (10X* 35HTHF MAXIMUM H E AT- T RAN S r E R RATE IS , 

1 El 2 * 5 » 7H WATTS// 

2 10X f 29HTHF TOTAL LIQUID IN NICK IS , 

3 £12.5,7H GRAMS/// 

4 l5X t 8HOISTANCE,6X,6HSTRtSS,7X f HHSTATIC HEA) t 6X* 9HPOROS IT X * 

5 5 X » 1 2 HCR IT. STRrE3S/17X,4H (CM* , 9X , 9H (CM LIG. > 9 6X, 9H (CM lIQ.I, 
b 22X,9H(GH LIQ.1//I 

966 FORMAT ( 1 1 X, 5E1 5 . 4). 

RETURN 

END 


tn 

VO 



